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SUMMARY 

A  modification  oj  the  doublet  lattice  method  of  A Ibano  and  Rodden  [Ref.  I]  has  been 
programmed  to  calculate  the  unsteady  generalised  airforces  and  pressures  acting  on  inter¬ 
fering  lifting  surfaces  in  subsonic  flow.  The  pre.sent  method  is  applicable  to  general  non- 
planar  and  nonparallel  lif  ting  .surfaces,  which  may  be  both  intersecting  and  nonintersecting 
in  nature.  The  relevant  theory  is  developed  in  detail  and  an  outline  of  the  program  and 
its  usage  is  given.  Comparisons  with  the  results  of  other  workers  are  included. 
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1.  INTRODUCTION 


The  aeroelastic  stability  of  an  aircraft  is  an  important  consideration  in  the  overall  process 
of  determining  its  airworthiness.  When  external  bodies  or  stores  are  added  to  the  wing  of  an 
aircraft,  the  dynamic  characteristics  of  the  aircraft  will  be  changed.  Due  to  the  inertial,  elastic 
and  aerodynamic  coupling  between  the  wing  and  its  stores  the  flutter  speed  of  the  aircraft  may 
be  adversely  affected.  The  aerodynamic  coupling  may  be  particularly  important  if  the  stores 
have  aerodynamic  surfaces,  such  as  fins,  that  can  generate  large  oscillatory  aerodynamic  forces. 

At  ARL  there  are  a  number  of  computer  programs  which  predict  unsteady  aerodynamic 
forces  acting  on  oscillating  lifting  surfaces.  These  programs,  which  are  based  on  the  doublet 
lattice  method  of  Albano  and  Rodden  [Ref.  1],  are  described  in  Reference  2  and  have  been 
optimised  for  simple  interfering  wing-tail-fin  and  control  surface  cases.  They  are  not  applicable 
to  general  combinations  of  interfering  lifting  surfaces. 

Thus  a  requirement  exists  for  a  computer  program  that  can  calculate  unsteady  aerodynamic 
forces  for  general  configurations  of  interfering  lifting  surfaces,  including  a  capability  for  dealing 
with  store  aerodynamics.  Since  the  doublet  lattice  method  is  readily  applicable  to  interfering 
and  nonplanar  lifting  surfaces,  it  is  used  as  the  basis  of  a  program  capable  of  dealing  with  general 
configurations.  The  doublet  lattice  method  is  versatile  and  there  are  essentially  no  restrictions 
on  the  configurations  that  can  be  handled  as  long  as  an  appropriate  idealisation  can  be  developed. 


2.  THE  DOUBLET  LATTICE  METHOD 

This  report  is  primarily  concerned  with  the  unsteady  aerodynamic  forces  generated  by  the 
oscillatory  motion  of  a  lifting  surface.  The  surface  is  idealised  as  a  thin  flat  plate  oscillating 
in  a  potential  flow.  The  steady  forces  associated  with  the  thickness  and  camber  of  the  lifting 
surface  are  not  considered. 

The  pressure  difference  which  exists  between  the  upper  and  lower  surfaces  of  the  plate  is 
integrated  to  give  the  aerodynamic  forces  of  interest.  This  oscillatory  pressure  is  calculated  by 
replacing  the  lifting  surface  by  a  planar  array  of  unsteady  pressure  doublets  which  have  strengths 
that  are  chosen  such  that  the  boundary  condition  of  tangential  flow  is  satisfied.  The  doublet 
strengths  are  related  to  the  induced  normalwash  by  an  integral  equation  which  for  general 
combinations  of  lifting  surfaces  must  be  solved  by  approximate  numerical  methods. 

In  Reference  1  Albano  and  Rodden  give  the  linearised  integral  equation  relating  the  induced 
oscillatory  normalwash  w{x,y,  to  the  pressure  distribution  p(x^,  yo,  over  all  lifting 

surfaces  L.S.  as: 


w'(.Y,  y,z)  1 
U 


K(x—Xfi,y  yo,z  —  zo',  M,  u))p(xo,  zo)dS 
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where  K  is  the  subsonic  nonplanar  Kernel  function,  which  is  a  function  only  of  geometry, 
Mach  number  M.  and  frequency  w.  The  symboljjc  indicates  integration  in  the  sense  of  Mangier 
[Ref.  3].  The  coordinates  of  the  sending  and  receiving  points  are  given  by  (xo,  yo,  co)  and 
(jY,  y.  z)  respectively. 

In  the  doublet  lattice  method  it  is  assumed  that  the  lifting  surface  can  be  approximated  by 
segments  of  planes.  The  surface  is  divided  into  small  trapczial  panels  with  parallel  sides  which  lie 
streamwise  and  leading  and  trailing  edges  which  are  some  linearly  varying  proportion  of  the 


local  chord  behind  the  leading  edge  of  the  surface  (Fig.  1).  The  panels  are  selected  such  that  any 
surface  discontinuities  such  as  control  surface  edges  or  surface  intersections  lie  on  panel  edges. 
Also,  the  panels  in  any  coplanar  streamwise  column  are  of  the  same  width. 

By  assuming  that  the  unknown  pressure  p  is  uniform  over  a  panel,  the  integral  of  Equation 
2.1  may  be  approximated  as  follows: 


u(.v,  z) 
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c  -  Cot  M,  iu)dS 
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/  Panel  J 


where  /  is  an  index  indicating  the  sending  panel.  The  integration  of  K  in  the  streamwise  direction 
is  done  simply  by  lumping  the  effect  into  a  loaded  line  of  pressure  doublets  at  the  1  /4-chord 
line  of  each  panel.  Hence  Equation  2.2  becomes: 


it(.v,  V,  r)  \  /Jj  (  jcos^j 
U  ~  2^  plf^  477 


^K{.\'  —  .Vo,  f  to.  -  ~co'.  M.  cl))JI 


vs  here 


/  /i 

(■]  is  the  mean  chord  of  the  /'th  panel, 

/3j  is  the  sweep  angle  of  the  1,4-chord  line  of  the  /'th  panel. 
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and  the  subscript  /j  denotes  that  the  line  integral  is  to  be  carried  out  over  the  1  /4-chord  line  of  the 
y'th  panel.  Note  that  for  any  panel  dS  =  c  cos  ^dl. 

In  Equation  2.3  the  normalwash  boundary  condition  ti(.v,  j.  r)  is  known  (see  Section  6) 
and  the  pressure  p,  over  each  panel  is  unknown.  A  set  of  linear  algebraic  equations  may  be 
formed  from  Equation  2.3  if  the  normalwash  is  satisfied  at  as  many  points  as  there  are  panels. 
There  is  one  control  or  receiving  point  per  panel  and  the  surface  normalwash  boundary  condition 
IS  satisfied  at  each  of  these  points.  The  control  point  is  located  at  the  mid-span  3 /4-chord  point 
of  each  panel. 


3.  MATRIX  FORM  OF  DOUBLET  LATTICE  METHOD 

As  shown  by  Albano  and  Rodden  [Ref  I),  the  doublet  lattice  method  reduces  the  integral 
equation  of  Equation  2.1  to  a  set  of  simultaneous  equations.  These  may  be  written  in  matrix 
form  as: 


[D]\y, 

where  is  a  complex  column  matrix.  The  /'th  clement  is, 

11  ,(.Vi.  IS.  Cl) 
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3.1 
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and  t,  is  the  nondimcnsional  normalwash  at  the  mid-span  3  4-chord  point,  (  vi.  i  ,.  c,). 
of  the  /'th  panel : 

[/)]  IS  a  complex  square  matrix  of  aerodynamic  influence  coefficients.  The  element  in 
the  /'th  row  and  /"th  column  is. 
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A, ,  dl 


and  Z)ij  is  the  line  integral  of  the  Kernel  function  along  the  1  /4-chord  line,  /j,  of  the 
y’th  panel ; 

{A}  is  a  complex  column  matrix.  The y’th  element  is, 


and  Aj  is  the  nondimensional  pressure  acting  uniformly  over  the  y’th  panel; 

Cj  is  the  mean  chord  of  the y’th  panel ; 

jSj  is  the  sweep  angle  of  the  1/4-chord  line  of  the  /'th  panel ; 

Kii  is  the  subsonic  nonplanar  Kernel  function  relating  the  normalwash  at  the  control  point 
on  the  /’th  panel  to  the  pressure  occurring  over  the  /’th  panel. 


4.  EVALUATION  OF  THE  KERNEL  FUNCTION 

4.1  Expression  for  the  Oscillatory  Kernel  Function 

In  Reference  4,  Landahl  has  shown  that  the  oscillatory  subsonic  nonplanar  Kernel  function 
K,  which  relates  the  complex  normalwash  h<x,  y,  to  the  complex  pressure  p(xo,  yo,  2o)e‘“*, 
can  be  written  in  the  following  form  : 

K  =  e U(A’i  Ti  +  K^nyr^  4. 1 

where 

Ti  =  cos  (yr  — ys) 

T2  =  {zi  cos  yr-yi  sin  yr}{zi  cos  y.s— >’i  sin  y„}/r^ 
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and  yr  and  y.s  are  the  dihedral  angles  of  the  receiving  and  sending  points  located  at  (  v,  r,  z) 
and  (jfo,  yo,  zo)  (see  Appendix  A  and  Figure  2).  The  terms  h  and  I2  represent  the  two  infinite 
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4.2  Evaluation  of  Ii  and  I2 


In  order  to  evaluate  the  Kernel  function  the  two  infinite  integrals  /i  and  />  must  be  deter¬ 
mined.  Since  they  cannot  be  expressed  in  terms  of  known  functions  they  must  be  evaluated 
numerically.  Computational  efficiency  is  a  prime  consideration  when  developing  a  procedure 
for  calculating  A  and  />  since  they  will  be  calculaud  numerous  times.  A  number  of  alternative 
procedures  for  evaluating  /i  and  A  will  now  be  described. 

In  Reference  1,  Albano  and  Rodden  integrate  Equation  4.2  by  parts  to  yield  : 


4.4 


The  integrand  in  the  last  term  of  the  above  equation  is  approximated  by  a  three-term  exponential 
function  and  the  integration  is  performed  analytically.  The  integral  A  is  evaluated  in  a  similar 
fashion,  requiring  integration  of  Equation  4.3  by  parts  twice. 

In  Reference  2,  Farrell  makes  use  of  the  expansions  given  by  Davies  [Ref.  5]  to  express 
A  and  /•>  as: 
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where  H  i  and  H  >  are  Struve  functions.  A  and  A  are  modified  Bessel  functions  of  the  first 
kind,  and  Ai  and  A>  are  modified  Bessel  functions  of  the  second  kind.  These  are  evaluated  by 
direct  series  summation  [Ref.  6].  with  the  asymptotic  expansions  being  used  for  larger  values 
of  r.  The  finite  integrals  of  Equations  4.5  and  4.6  are  evaluated  by  a  numerical  technique  which 
divides  the  range  of  integration  into  a  number  of  regions,  the  lengths  of  which  are  functions 
of  i>.  The  integration  is  performed  over  each  of  these  smaller  regions  using  six-point  Gaussian 
quadrature.  For  values  of  iii  that  are  large  in  magnitude  the  expressions  for  /i  and  A  simplify 
considerably  since  many  of  the  terms  may  be  neglected.  Although  this  method  is  more  accurate 
than  that  of  Reference  I.  the  computation  time  is  con,siderably  greater. 

In  Reference  7.  Geising  cf  at.  improve  on  the  method  used  in  Reference  I  by  replacing 
the  three-term  exponential  approximation  to  the  integrand  of  Equation  4.4  with  an  eleven-term 
exponential  approximation  due  to  Laschka  [Ref.  8).  The  maximum  error  of  this  approximation 
is  O'  I3.‘'"„  [Ref.  7).  and.  as  before,  its  use  allows  l\  and  A.-  to  be  evaluated  analytically.  A  com¬ 
parison  between  this  method  ;  nd  that  of  Reference  2  was  carried  out  by  calculating  trial  values 
of  the  Kernel  function.  Good  agreement  between  the  two  approaches  was  noted,  and  the  method 
of  Reference  7  was  adopted  since  it  reduced  the  computation  time  involved  in  the  calculation 
of  the  Kernel  function,  as  well  as  yielding  results  of  a  high  accuracy. 

A  detailed  description  of  the  method  used  for  calculating  h  and  A  is  to  be  found  in  Appen¬ 
dix  B, 
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4.3  Expression  for  the  Steady  Kernel  Function 

The  expression  for  the  steady  Kernel  function  given  in  Reference  9  nray  be  written  in  the 
following  form : 

Ko  =  (/ri'S'n  +  K2f^'T2)/r^  4.7 

where 

Ti  =  cos  (yr—ye) 

Tz  =  [zi  cos  yr-y\  sin  cos  ys->’i  sin  ys}/r2 


and  the  other  symbols  are  as  defined  in  Equation  4.1.  The  expression  for  the  steady  Kernel  is 
considerably  simpler  than  that  for  the  oscillatory  Kernel  function.  Hence  when  the  frequency 
of  oscillation  can  be  assumed  to  be  zero  it  is  advantageous  to  use  Equation  4.7  in  preference  to 
Equation  4. 1 . 


5.  INTEGRATION  OF  THE  KERNEL  FUNCTION 

At  this  point  it  is  desirable  to  discuss  the  nonplanar  Kernel  function  as  contrasted  with 
the  planar  one,  since  the  nature  of  the  Kernel  function  will  have  a  direct  bearing  on  the  integration 
scheme  that  is  adopted  (see  Appendices  C  and  D). 


5.1  Behaviour  of  the  Steady  Kernel  Function 

It  is  convenient  to  confine  our  investigation  to  the  behaviour  of  the  steady  Kernel  function 
since  this  function  is  amenable  to  an  analytical  approach.  The  underlying  assumption  of  this 
analysis  is  that  the  salient  features  of  the  steady  Kernel  function  are  representative  of  the  behaviour 
displayed  by  the  oscillatory  Kernel  function. 

Following  the  approach  of  Berman  et  al.  [Ref  10],  Figure  3  provides  a  comparison  of  the 
steady  nonplanar  Kernel  function  (Equation  4.7  with  zi  /  0)  and  the  steady  planar  Kernel 
function  (Equation  4.7  with  z\  =  0).  In  Figure  3  the  values  of  the  Kernel  function  are  plotted 
against  the  spanwise  variable  yi  in  the  vicinity  of  y\  =*  0  for  several  values  of  vertical  distance  zi 
and  for  chordwise  distance  .vi  1.  The  chordwise  variation  is  not  shown  since  the  nonplanar 
Kernel  function  is  very  similar  to  the  planar  Kernel  function  for  the  .xi  coordinate  [Ref  10). 

The  figure  clearly  shows  that  the  spanwise  variation  of  the  nonplanar  Kernel  function  in 
the  vicinity  of  the  control-point  station  is  very  different  from  that  of  the  planar  Kernel  function. 
The  singularity  at  yi  =  0  is  no  longer  present.  Instead,  a  sharp  reversal  in  the  trend  of  the 
Kernel  function  is  experienced,  resulting  in  a  large  negative,  but  fini'e,  value  at  yi  =  0.  Thus 
the  Mangier  concept  [Ref.  3]  of  the  “finite  part  of  infinite  integrals"  (cf  Section  2)  is  no  longer 
required  for  the  integration  of  the  nonplanar  Kernel  function. 
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In  Reference  10,  Berman  et  al.  found  that  the  twin  positive  peaks  of  the  nonplanar  Kernel 
function  are  located  at  >>1  =  ±  ziy/3  where  zi  is  the  vertical  separation  between  the  nonplanar 
(but  parallel)  lifting  surfaces.  Since  the  large  gradients  of  this  Kernel  function  occur  between 
these  peaks,  it  is  essential  that  the  relatively  small  region  ~zi\/3<}’i<  +  zi\/3  be  provided 
with  a  high  density  of  integration  stations. 

Although  the  function  shown  in  Figure  3  is  symmetric  about  >i  =  0,  this  is  not  the  case 
when  dealing  with  general  combinations  of  nonplanar  nonparallel  lifting  surfaces.  The  amount 
of  asymmetry  which  arises  is,  however,  small  in  comparison  with  the  symmetric  part  of  the 
Kernel. 

The  above  mentioned  behaviour  must  be  adequately  accounted  for  in  the  line  integral  of 
Equation  3.2.  If  suitable  integration  schemes  are  not  implemented  then  the  aerodynamic  in¬ 
fluence  coefficients  will  lose  accuracy  as  the  vertical  separation  between  any  two  surfaces  is  re¬ 
duced.  Finally,  at  very  small  values  of  vertical  separation,  as  ri->0,  the  results  will  be  completely 
inaccurate. 

The  previous  discussion  applies  to  nonplanar  horizontal  lifting  surfaces.  For  the  general 
case  of  combinations  of  nonplanar  nonparallel  lifting  surfaces  it  is  necessary  to  consider  the 
variable  rather  than  the  variable  zi  in  isolation.  For  this  general  case  numerical 

difficulties  in  the  integration  will  arise  when  r'^->0. 


5.2  Behaviour  of  Kor-  for  Nonplanar  Panel  Combinations 


it  is  clear  that  the  ill-conditioning  of  the  integrand  in  Equation  3.2  presents  computational 
difficulties.  By  analytically  incorporating  the  1  /r-  behaviour  of  the  Kernel  function  into  the  line 
integral  we  may  write  Equation  3.3  in  the  following  form  ; 


Cj  cos 


This  has  the  effect  of  minimising  the  variations  that  must  be  dealt  with  by  the  integration  scheme. 

In  order  to  maximise  the  accuracy  of  integration  it  is  worthwhile  examining  the  behaviour 
of  the  numerator  of  Equation  5.1  for  commonly  occurring  situations.  Once  again  we  confine 
our  attention  to  the  steady  Kernel  function  by  utilising  the  assumptions  presented  in  Section  5.1. 
For  the  steady  case  the  numerator  of  Equation  5.1  (neglecting  (  j  cos/3j/4w)  may  be  written 


k 


Kof^  =  + 


5.2 


The  typical  behaviour  of  Kor-  is  illustrated  in  Figure  4.  Values  of  Kor'^  are  plotted  against  the 
spanwise  variable  vi  in  the  vicinity  of  vi  =  0  for  three  values  of  vertical  distance  zi  and  a 
chordwise  distance  .vi  =  I.  with  both  the  sending  and  receiving  panels  lying  in  planes  parallel  to 
the  horizontal  plane.  Note  that  the  variation  in  the  vicinity  of  vi  =  0  is  now  considerably  less  than 
that  of  the  corresponding  curves  in  Figure  3.  Furthermore,  the  curves  are  asymptotic  to  the 
value  of  +10  for  values  of  ,vi-*  +  oo,  rather  than  zero  as  in  Figure  3. 

The  occurrence  of  aerodynamic  surfaces  that  are  perpendicular  to  each  other  is  quite  com¬ 
mon  (e.g.  a  T-tail  and  the  cruciform  fins  at  the  rear  of  a  missile).  For  such  cases  the  variation 
of  Kor-  differs  considerably  from  that  illustrated  in  Figure  4.  The  behaviour  of  KoF'^.for  perpen¬ 
dicular  surfaces  is  illustrated  in  Figure  5.  It  is  seen  that  there  is  only  one  zero  crossing  in  contrast 
to  the  two  zero  crossings  evident  in  Figure  4.  Also,  there  are  only  two  turning  points  instead 
of  three. 


5.3  Prediction  of  Zero  Crossings  and  Turning  Points 

Although  the  numerator  of  the  integrand  in  Equation  5.1  is  much  better  behaved  than  the 
integrand  of  Equation  3.2,  the  highly  localised  nature  of  Koi"'^  is  still  evident  in  Figures  4  and  5. 
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In  order  to  maximise  the  accuracy  of  integration  it  is  essential  that  the  relatively  small  region  of 
rapid  variation  be  provided  with  a  high  density  of  integration  stations.  Hence  it  is  necessary 
to  develop  procedures  for  determining  the  location  of  this  small  region  and  deciding  whether 
it  lies  on  the  interval  of  integration  /j  (defined  in  Equation  3.3). 

Following  the  approach  of  Berman  [Ref.  10]  an  attempt  was  made  to  derive  an  analytical 
technique  for  determining  the  turning  points  of  Kor^,  since  these  define  the  region  of  interest. 
This  particular  approach  was  abandoned  because  of  its  complexity.  Also,  the  long  compu¬ 
tation  times  would  have  made  any  calculations  relatively  expensive. 

A  considerably  simpler  and  more  efficient  method  resulted  when  attention  was  focussed 
on  the  zero  crossings  of  Kor'^,  rather  than  its  turning  points.  It  is  possible  to  write  yi  and  zi 
in  terms  of  a  parameter  /,  where  the  range  0  <  /  <  1  represents  the  interval  of  integration  /j. 
Hence,  for  any  such  interval,  the  Kernel  function  may  be  defined  as  a  function  of  /.  Turning 
our  attention  to  Figure  6,  the  left  and  right  zero  crossings,  /l  and  /r,  are  estimated  using  the 
method  described  in  Appendix  E.  The  turning  point,  /t,  is  then  located  by  numerical  sub¬ 
division  of  the  interval  between  the  zero  crossings.  Hence  the  two  intervals,  A/l  and  A/r,  may 
be  calculated  from ; 


A/l  =  /t— /l 

5.3 

A/r  =  /r— /t 

5.4 

Note  that  Figure  4  illustrates  the  special  case  where  AtL  =  A/r  and  /t  =  0-5,  since  in  general 
A/l  A/r  and  /t  does  not  occur  in  the  centre  of  the  interval  /j.  From  a  further  inspection  of 
Figure  4  it  would  appear  that  the  region  of  greatest  variation  of  Kor^  lies  approximately  within 
the  interval : 

/t — 6A/l  <  /  <  /t-1-6A/r  5.5 

Hence  this  region  should  be  provided  with  a  high  density  of  integration  stations. 

The  above  method  is  not  valid  for  aerodynamic  surfaces  that  are  perpendicular  to  each 
other  since  only  one  zero  crossing  exists  for  such  cases.  Due  to  a  number  of  simplifications  that 
occur,  it  is  possible  to  estimate  the  positions  of  the  turning  points  for  cases  where  the  surfaces 
are  perpendicular.  Referring  to  Figure  7,  the  left  and  right  turning  points,  /tl  and  /tr,  are 
estimated  using  the  method  described  in  Appendix  F,  whereas  the  zero  crossing,  to,  may  be 
estimated  using  the  method  of  Appendix  E.  Hence  we  may  define  the  two  left  and  right  intervals, 
A/tl  and  A/tr,  using : 

A/tl  =  to—tTL  5.6 

A/tr  =  /tr— to  5.7 

Clearly,  the  region  of  rapid  variation  falls  within  the  interval : 

/tl  <  to  <  /tk  5.8 

Once  again,  this  region  should  be  provided  with  a  high  density  of  integration  stations. 


6.  NORMALWASH  AND  PRESSURE  DISTRIBUTIONS 

In  setting  up  a  system  of  panels  to  model  an  actual  case  it  is  assumed  that  the  system  is 
capable  of  vibration  in  a  number  of  modes  of  displacement.  In  mode  p.  the  displacement  normal  to 
the  mean  position  of  a  point  (.v,  y,  r)  on  the  surface  is  taken  to  be  / /p(,v,  v,  z).  When  the  system  is 


oscillating  harmonically  about  its  mean  position,  with  circular  frequency  w,  in  the  mode  p, 
the  displacement  normal  to  the  surface  is  taken  to  be 


Zp(.v,  >■,  r,  t)  =  y, 


6.1 


where  /  is  the  reference  length. 

Applying  the  boundary  condition  that  the  flow  must  be  tangential  to  the  surface  requires 
the  relationship  between  the  normalwash  and  the  modal  displacement  to  be 


u  p(.v,  ,1’. 


DZp 

Dt 
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This  yields  the  following  expression  for  the  nondimensional  normalwash 


up(.Y,  V,  :) 

U 


a(.v/7) 


/p(.v,  j',  z)  +  ik/p(.\\y,  :) 


=  a'  +  ik<x" 


6.3 

6.4 


where  k  is  the  frequency  parameter  cul/U.  Thus  the  real  part  a'  of  the  normalwash  matrix  {«} 
of  Equation  4.1  contains  the  chordwise  slopes  of  the  displacement  mode  at  the  control  point 
in  each  panel  and  the  scaled  imaginary  part  a"  contains  the  modal  displacements  at  these  points. 

When  the  normalwash  influence  coefficient  matrix  [D]  has  been  calculated  and  the  non- 
dimensional  normalwash  {a|  has  been  prescribed,  the  system  of  linear  equations, 

{«}  -  [0]{A}, 


defined  previously  in  Equation  3.1,  may  be  solved  to  find  the  nondimensional  pressure  {A}  acting 
over  all  surfaces.  Note  that  there  are  as  many  lifting  pressure  distributions  as  there  are  modes. 


7.  GENERALISED  FORCES 


For  dynamical  analyses  of  the  vibration  of  aerodynamic  surface  combinations  we  generally 
apply  Lagrange's  equations  of  motion.  This  requires  expressions  for  the  generalised  airforces 
that  occur.  These  airforces  act  on  the  actual  aerodynamic  surfaces  but  within  the  linearised 
approximation  we  can  take  the  pressures  described  above  and  evaluate  the  generalised  airforces. 
The  generalised  force  Fp,,  in  the  mode  p  due  to  oscillation  in  the  mode  q  is  given  by 


P  p(| 


r  r 


//p(,v,  .1’,  z)pU-A,,(.x.  y.  z)ds  e'' 


L.S. 


pUV  M-V.  v,  r)A„(.v,y,  z)(iS 
L.S. 

pU-r-^Qp.,e"'"  7.1 

where  L.S.  represents  the  integral  over  all  the  lifting  surfaces,  /p  is  the//th  nondimensional  dis¬ 
placement  mode.  A,,  is  the  nondimensional  pressure  due  to  the  r/’th  mode  and  k  is  the  frequency 

lu/ 

parameter 


From  Equation  7. 1  the  nondimensional  generalised  force  Qpq  is 


Qm  - 


fxix,  y,  z)Aq(jf,  y,  z)dS 


L.S. 


1.2 


In  the  doublet  lattice  method,  where  the  aerodynamic  surfaces  are  approximated  by  panels 
and  the  pressure  is  assumed  to  be  constant  over  each  panel,  Equation  7.2  simplifies  to 


n 


Qp<\  — 


z 


;•  =  1 


7.3 


where  /p‘  is  the  nondimensional  displacement  at  the  mid-span  quarter-chord  point  of  the  /’th 
panel  in  the  ^j’th  mode,  Aq‘  is  the  pressure  on  the  /'th  panel  due  to  the  ^'th  mode,  l-A  -,  is  the  area  of 
the/’th  panel  and  n  is  the  total  number  of  panels  used  to  model  the  aerodynamic  surfaces.  The  mid¬ 
span  quarter-chord  point  is  called  the  lift  point  and  is  illustrated  in  Figure  I.  Finally,  we  note 
that  it  is  customary,  for  dynamical  analyses,  to  write  in  the  form 

Qp<\  —  Q'p(\  +  il<Q"p(i  7.4 

where  0'pq  and  5  no  ^re  real  quantities. 


8.  IDEALISATION  OF  GENERAL  CONFIGURATIONS 

In  order  to  minimise  the  amount  of  time  spent  in  developing  a  panel  distribution  to  idealise 
a  particular  configuration,  it  is  essential  to  use  a  procedure  that  defines  the  panel  distribution  by 
interpreting  a  concise  and  compact  data  set.  Clearly,  it  is  not  desirable  to  h^'  e  to  define  each 
panel  individually  as  this  would  be  very  time  consuming  and  the  likelihooi  if  errors  would 
be  high. 


8.1  Definition  of  Lifting  Surface  Groups 

A  lifting  surface  group  is  composed  of  one  or  more  panels  and  it  is  a  subset  or  element  of 
the  idealised  configuration.  For  our  purposes  it  is  looked  upon  as  being  the  basic  unit  or  building 
block  which  allows  us  to  generate  sequences  of  panels  compatible  with  the  requirements  of  the 
doublet  lattice  method. 

A  group  is  defined  by  reference  to  four  nodes.  Figure  8  illustrates  the  node  numbering 
system  associated  with  the  topology  of  a  typical  group.  Nodes  1  and  2  define  the  leading  edge 
of  the  group,  and  nodes  3  and  4  specify  the  trailing  edge.  This  node  numbering  convention 
defines  the  topology  of  the  group  and  must  be  strictly  adhered  to.  Furthermore,  nodes  I  and  3 
must  lie  on  the  same  streamwise  line.  This  requirement  must  also  be  satisfied  by  nodes  2  and  4. 

It  is  useful  to  transform  the  trapezial  group  onto  a  square  defined  in  a  nondimensional 
coordinate  system,  wherein  the  coordinates  of  nodes  I.  2,  3  and  4  are  as  illustrated  in  Figure  9. 
Node  I  is  taken  as  the  origin  of  the  nondimensional  coordinate  system  associated  with  each 
group.  The  chordwise  and  spanwisc  ordinates,  i,-  and  t-,  are  positive  in  the  directions  shown. 
They  have  been  nondimensionaliscd  with  respect  to  the  local  chord  and  total  span  of  the  group. 
Note  that  AC  and  A5  refer  to  the  number  of  chordw  ise  rows  and  spanwise  columns  of  panels 
pertaining  to  the  group.  To  facilitate  panel  generation,  it  is  necessary  to  provide  only  the  non- 
dimensional  chordwise  and  spanwise  ordinates  of  the  edges  of  panels  in  a  particular  group. 
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The  panels  within  a  given  group  are  numbered  1,  2,  3,  .  .  .,  NC*NS  in  a  columnwise  pro¬ 
gression.  As  shown  in  Figure  9,  the  number  of  the  panel  in  the  /’th  chordwise  row  and  J'ih 
spanwise  column  is  given  by  the  following  equation : 

panel  number  =  (7— l)*jVC+ /  l^/^AC  8.1 

Note  that  the  chordwise  rows  of  panels  do  not  have  to  be  a  constant  proportion  of  the  local 
chord  across  the  span  of  the  group ;  a  linear  variation  is  acceptable.  Hence  it  is  possible  to  have 
one  or  more  rows  of  panels  of  constant  (dimensional)  chord  defined  within  a  group  that  is  tapered. 
Some  of  the  possible  combinations  are  illustrated  in  Figure  10.  It  is  also  possible  to  define 
evenly  spaced  panels  in  either  the  chordwise  or  spanwise  directions,  or  both,  without  having  to 
specify  the  ordinates  along  the  sides  of  the  group. 


8.2  Pressure  and  Normalwash 

The  definition  of  the  direction  of  positive  pressure  and  positive  normalwash  acting  on  a  set  of 
panels  in  a  group  is  implicit  in  the  topology  of  the  group.  Since  all  panels  in  a  group  share  the 
topology  of  that  group,  the  sense  of  positive  pressure  and  normalwash  is  likewise  defined.  Note 
that  the  normalwash  is  positive  in  the  same  sense  as  the  topology-defined  positive  pressure. 
Furthermore,  the  sense  of  positive  pressure  can  be  reversed  from  that  defined  by  the  topology  of 
the  group. 

Figure  1 1  illustrates  the  pressure  p  and  normalwash  ti-  acting  on  a  panel.  The  pressure  and 
normalwash  are  shown  acting  in  their  positive  sense.  As  mentioned  above,  the  topology  of  the 
panel,  represented  by  the  node  numbers  I,  2,  3  and  4.  echoes  the  topology  of  the  group  to  which 
It  belongs. 

Figure  12  shows  a  view,  looking  forward,  of  some  possible  orientations  of  panels.  The 
numbers  1  and  2  refer  to  the  topology  of  the  group  from  which  the  panels  derive.  The  arrows 
normal  to  the  surface  of  each  panel  indicate  the  sense  of  positive  pressure  as  defined  by  the 
topology  of  the  panel.  Although  no  vertical  or  horizontal  panels  are  shown,  their  characteristics 
are  readily  deduced  by  simply  projecting  the  existing  panels  onto  the  desired  axis. 


8.3  .Symmetry  and  Reflection  Planes 

The  majority  of  applications  of  the  doublet  lattice  method  involve  configurations  with  one 
or  more  planes  of  symmetry.  The  X-Z  plane  is  a  plane  of  symmetry  for  most  cases  ;  the  right 
side  of  a  lifting  surface  is  the  mirror  image  of  the  left  side.  In  other  cases  an  additional  plane  of 
symmetry  (the  X~  Y  plane)  exists.  For  instance  a  lifting  surface  in  the  proximity  of  the  ground 
(aircraft  in  ground  effect)  represents  a  configuration  with  two  planes  of  symmetry. 

In  general,  all  flow  conditions  can  be  split  up  into  symmetrical  and/or  antisymmetrical 
parts  in  relation  to  cither  plane  of  symmetry.  If  the  configuration  is  symmetrical  and  the  flow 
is  either  symmetrical  or  antisymmetrical  then  considerable  savings  in  computational  effort  can 
be  realised.  These  are  apparent  as  a  reduction  in  the  number  of  normalwash  influence  coefficients 
that  need  to  he  calculated  and  stored  in  computer  memory.  This  minimises  the  size  of  the  system 
of  equations  that  needs  to  be  solved. 

In  order  to  make  use  of  any  possible  simplifications  due  to  symmetry  or  antisymmetry  it  is 
necessary  to  he  able  to  create  a  system  of  image  lifting  surfaces  by  reflecting  groups  in  the  -V  Z 
and  \  )  planes.  Consider  the  general  case  where  a  lifting  surface  configuration  is  made  up  of  a 
source  group  and  three  additional  images.  This  is  illustrated  in  Figure  I  .C  Once  again  the 
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numbers  1  and  2  refer  to  the  topology  of  the  group  and  the  arrows  indicate  the  sense  of  positive 
pressure.  The  superscripts  (1),  (2),  (3)  and  (4)  identify  how  the  group  was  obtained  as  follows : 


(1)  — source  group  from  which  all  images  are  developed  ; 

(2) — image  obtained  by  reflecting  the  source  group  in  the  X-Z  plane ; 

(3)  — image  obtained  by  reflecting  the  source  group  in  the  X-Z  plane  followed  by  a  reflection 

in  the  X-  Y  plane ; 

(4)  — image  obtained  by  reflecting  the  source  group  in  the  X-  Y  plane. 

Note  that  the  topology  of  the  image  groups  in  relation  to  the  source  group  follows  the  con¬ 
vention  established  in  Figure  12.  Care  must  be  exercised  in  using  the  image  capability  as  anomo- 
lous  configurations  can  arise  under  some  conditions. 

Considering  the  most  general  case  of  a  source  and  three  image  groups,  as  depicted  in  Figure 
13,  the  matrix  equation  given  by  Equation  3.1  can  be  written  as : 


vv'i’  \ 

dll  di2  di3  di4 

i  ] 

H.12)  I 

dzi  dzz  doa  da^ 

1 

W'3>  1 

dai  daa  daa  da4 

L  ^^41  d42  d^a  dm 

1  ) 

where  the  i/ij  are  submatrices  of  matrix  D. 

For  symmetric  and  antisymmetric  cases  the  magnitudes  of  /t'®’  and  p*'*’  are  equal 

and  only  their  phases  can  differ  by  180"’  relative  to  say.  Hence  we  may  write : 


'  pi>) 
p(2l 

p'3’ 

pi4, 


S'" 

8<2I 

§(3I 

8<4) 


8.3 


where  the  terms  S'"  may  take  the  following  values 

5<i.  =  +1,  _| 

5<2)  =  +1,  _|,o 

5'3'  =  +1,  1.0 

§14)  =4-1,  -1,0 

Since  it  refers  to  the  source  group,  0,  whereas  since  some  image  groups  may  not  be 

included  in  a  given  configuration  the  appropriate  S  terms  will  be  zero.  If  8'"  =  1  then  the 
pressure  is  positive  in  the  same  sense  as  that  defined  by  the  topology  of  the  source  or  image  group. 
If  8*"  ^  1  then  the  pressure  is  positive  in  the  opposite  sense  to  that  defined  by  the  topology. 

By  substituting  Equation  8. .3  into  Equation  8.2  we  obtain  : 

m"’--  [8"V/n4-8'-’'r/|.4-8'3'f/i3+8'''>r/i,]/j'»  8.4 

Hence  the  pressure  distributions  for  all  of  the  groups  have  been  reduced  to  one  distribution  for 
the  source  group.  Thus  only  that  part  of  the  interaction  that  involves  the  source  group  needs 
to  be  considered,  since  the  others  may  be  readily  deduced  from  the  symmetry  or  antisymmetry 
conditions. 

Whether  it  be  symmetric,  antisymmetric  or  asymmetric,  each  type  of  case  is  defined  explicitly 
by  the  chosen  orientation  of  the  positive  senses  of  pressures  for  the  source  and  image  groups. 

1 1 


T* 


pvwwn 


MU^V«U«V« 


'■j  «u  «  u  V  u  «u  V 


The  topology-defined  senses  are  altered  by  the  appropriate  choice  of  values  for  the  8  factors  in 
Equation  8.4.  An  alternative  interpretation  is  to  view  S'>*  as  an  overall  multiplier  determining 
the  senses  of  the  reference  source  pressures  The  remaining  factors  S'2>,  8(3)  and  S‘‘*'  then 
determine  whether  the  effects  of  the  image  surfaces  are  to  be  added  to  or  subtracted  from  the 
source  normal  wash  factor  dn  in  order  to  obtain  the  correct  realisation  of  symmetry  or  anti¬ 
symmetry. 
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8.4  Example  of  Symmetry  and  Reflection  Planes 

As  an  example  of  the  concepts  described  in  Section  8.3,  let  us  idealise  a  model  of  a  T-tail 
in  a  wind  tunnel.  It  is  assumed  that  the  fin  and  the  floor  of  the  wind  tunnel  provide  two  axes  of 
symmetry.  The  effects  of  the  remaining  wind  tunnel  boundaries  are  neglected.  For  symmetric 
and  antisymmetric  modes  only  two  source  groups  need  be  defined,  the  effects  of  the  other  surfaces 
being  obtained  by  the  use  of  images.  Figure  14  shows  the  idealisation  of  the  T-tail  where  it  is 
assumed  that  the  senses  of  pressures  /Ji'*’  and  /)■>"’  define  the  topology  of  the  source  groups. 

In  order  to  model  symmetric  motion  of  the  T-tail  we  have  that  8)'-'  —  - 1,  Si'^’  -  -t- 1, 
8)'^'  ~  -I  (assuming  that  5i'"  =  -i-l).  Since  the  fin  does  not  sustain  a  pressure  difference 
8o'>'  and  8>'^'  are  both  equal  to  zero,  which  means  that  the  fin  and  its  image  do  not  need  to  be 
included  in  the  analysis.  Keeping  in  mind  the  fact  that  ^2"’  =  0,  Equation  8.4  may  be  written  as ; 


(  "'i'*'  I  _  ( +  I  Wu -l-(  —  l)tfi2-(-( -h  l)^ft3  +  (  —  l)^fl4  (  pi'*’  1 

I  U'2**'  /  (0)d21+  (0)<f22+  (0)^f23+  (0)fl'24  (  ^2*”  j 


8.5 


where  the  fact  that  82“’  =  0  indicates  that  the  problem  may  be  simplified  by  omitting  the  fin 
from  the  symmetric  analysis. 

For  the  antisymmetric  case  Equation  8.4  becomes: 


(-(- l)r/n  +  (+ l)<fi2-t-(  — l)<fi3  +  (  — l^u  I  pi'” 

(+I)d2l+  (0)d22+  {0)d23  +  (~  l)d2A  | 


8.6 


Hence  it  is  seen  that  symmetric  and  antisymmetric  cases  can  be  dealt  with  by  appropriately 
specifying  the  applicable  8  factors. 


9.  COMPARISON  WITH  OTHER  METHODS 

In  order  to  illustrate  the  method  di.scu.ssed  above,  which  has  been  implemented  at  ARE 
into  the  FORTRAN  computer  program  AIRFORTP,  three  particular  interfering  lifting  surface 
configurations  have  been  investigated  over  a  range  of  frequency  parameters,  Mach  numbers 
and  variations  in  geometry.  These  configurations  are  described  in  Sections  9.1  to  9.3  and  the 
results  compared  vvith  those  obtained  using  alternative  approaches  are  shown  in  Tables  I  to  9. 
The  primary  interest  is  the  comparison  of  the  results  from  the  various  computational  formula¬ 
tions  of  the  doublet  lattice  method  and  the  lifting  surface  results  are  included  only  as  a  secondary 
comparison. 

Note  that  in  all  the  Tables  the  generalised  forces  of  the  form  0,,,,  are  presented  in  terms  of 
modulus  and  phase  The  phase  angle  is  given  in  degrees  and  lies  in  the  range 

0  ^  ^  ,160  . 
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9.1  AGARD  Horizontal  Wing  and  Tailpiane 


The  AGARD  horizontal  wing  and  tail  combination  is  shown  in  Figure  15.  This  is  the 
configuration  which  AGARD  specified  for  calculation  of  generalised  airforces  and  the  wing  and 
tailpiane  are  swept-back  tapered  wings  in  close  proximity.  The  semi-spans  of  the  wing  and  tail- 
plane  are  both  of  unit  length  and  the  reference  length  is  taken  equal  to  the  semi-span.  The  origin 
of  coordinates  is  at  the  apex  of  the  wing. 

The  configuration  is  in  a  subsonic  flow  of  free-stream  Mach  number  M  =  0  •  8  and  is 
oscillating  with  a  frequency  parameter  k  in  the  range  0  to  1  -  5  in  one  of  two  antisymmetric  modes 
of  oscillation,  defined  analytically  by 


f\(x,y,  z)  =  >’(Ar-2-25  |3’ !  -0-85) 

on  the  wing 

9.1 

=  >' 

on  the  tail 

9.2 

Hx,  y,  z)  =  y\  y  ! 

on  the  wing 

9.3 

=  (Af-3-35)sg«(>) 

on  the  tail 

9.4 

The  first  mode  is  torsion  of  the  wing  about  an  axis  38%  of  the  local  chord,  coupled  with  roll  of 
the  tail.  The  second  is  parabolic  bending  of  the  wing  coupled  with  pitch  of  the  tail. 

The  configuration  with  /;  #  0  is  obtained  from  the  configuration  when  the  wing  and  tail- 
plane  are  coplanar  simply  by  translating  the  tailpiane  a  distance  h  parallel  to  the  z-axis. 

The  comparison  of  results  obtained  for  the  AGARD  configuration  is  largely  based  on  data 
collected  and  presented  by  Davies  [Ref.  5].  Tables  1  to  4  compare  the  results  published  in 
Reference  5  with  the  results  of  Farrell  [Ref.  2]  and  the  present  theory.  The  results  attributed 
to  Davies  and  Albano,  Perkinson  and  Rodden  [Ref.  II]  are  based  on  lifting  surface  methods, 
the  remainder  being  based  on  variations  of  the  doublet  lattice  method. 

In  the  application  of  the  present  method  an  array  of  8  panels  along  the  semi-span  and  8 
panels  along  the  chord  was  used  for  both  the  wing  and  the  tailpiane.  These  panel  distributions 
were  evenly  spaced  in  the  spanwise  and  chordwise  directions  As  a  result  the  panels  on  the  wing 
and  the  tail  were  aligned  in  streamwise  strips. 

Table  I  presents  generalised  airforces  for  two  values  of  frequency  parameter,  ^  =  0  and 
k  =  1  5,  for  the  case  where  the  wing  and  tailpiane  are  coplanar  (/;  =  0).  Inspection  of  Table  I 
shows  that  the  present  method  compares  favourably  with  the  results  of  the  other  doublet  lattice 
and  lifting  surface  methods.  Particularly  good  agreement  is  obtained  with  the  refined  doublet 
lattice  method  of  Kalman,  Geising  and  Rodden  [Ref.  12).  However,  for  /c  =  1  -  5  the  present 
method  consistently  over-estimates  the  magnitudes  of  the  generalised  forces  by  an  average 
of  10  •6'’,,  when  compared  with  the  results  obtained  by  Farrell,  although  the  calculated  phase 
angles  are  in  good  agreement.  This  is  in  contrast  to  the  generally  good  agreement  obtained 
for  k  -  0. 

Table  2  presents  resul's  for  generalised  airforces  for  the  nonplanar  case  with  h  =  0-6  for 
frequency  parameters  ^  =  0  and  A'  —  1-5.  Once  again  the  results  for  the  present  method  agree 
favourably  with  those  of  other  workers.  Previously  made  comments  regarding  the  comparison 
of  the  present  method  and  that  due  to  Farrell  for  the  case  where  A  ^  I  -  5  are  equally  applicable 
to  the  results  in  Table  2. 

Table  3  compares  Davies  lifting  surface  results  for  A  0  6  and  a  range  of  frequency  para¬ 
meters  from  A  =  0  to  A  -  15  with  results  obtained  by  the  present  method.  The  main  purpose 
of  the  comparison  is  to  check  the  variation  of  the  generalised  forces  with  frequency  parameter. 
The  agreement  between  the  trends  predicted  by  the  two  methods  is  satisfactory.  The  variation 
between  the  results  in  this  comparison  is  within  the  variation  exhibited  by  the  results  shown  in 
the  previous  Tables. 

Table  4  compares  Davies  lifting  surface  results  for  A  1-5  and  a  range  of  vertical  separa- 
tionsfrom/i  OtoA  0-6  with  results  obtained  by  the  present  method.  The  results  for /i  =  001 
and  h  0  04  are  of  particular  interest  since  they  represent  a  nonplanar  case  with  a  very  small 
vertical  separation  between  the  wing  and  tailpiane.  Small  vertical  separations  between  surfaces 
can  lead  to  problems  with  the  doublet  lattice  method  such  as  those  discussed  in  Section  5.1. 


Comparison  of  the  present  method  with  that  of  Davies  indicates  that  there  is  good  agreement 
in  the  trends  predicted  for  the  generalised  airforces  as  the  vertical  separation  is  varied,  together 
with  an  acceptable  match  between  individual  pairs  of  generalised  forces.  In  the  critical  region 
of  interest  when  li  is  less  than  01  there  is  acceptable  agreement  between  the  two  methods. 


! 

‘  9.2  Stark’s  Swept  and  Tapered  T-Tail 

I  A  swept  and  tapered  T-tail,  first  analysed  by  Stark  in  Reference  13.  is  shown  in  Figure  16. 

'  A  feature  of  this  particular  T-tail  configuration  is  that  at  the  stabiliser-fin  junction  the  trailing 

I  edge  of  the  stabiliser  extends  beyond  the  trailing  edge  of  the  fin.  The  reference  length  for  this 

I  configuration  was  taken  to  be  the  stabiliser  semi-span. 

j  The  T-tail  is  oscillating  in  three  rigid  body  modes.  These  consist  of  yawing  about  a  vertical 

,  axis  through  the  centre  of  the  root  chord  of  the  fin  (positive  nose  right),  sidesway  (positive  left) 

'  and  rolling  about  the  fin-stabiliser  intersection  (positive  right  stabiliser  down).  These  modes 

,  may  be  expressed  analytically  as 


/|(.V.  ,r,  r) 

.  3(.v-i-0- 15577) 

on  the  fin 

9.5 

-  0 

on  the  stabiliser 

9.6 

/L'(.v.  .1'.  --) 

.  10 

on  the  fin 

9.7 

0 

on  the  stabiliser 

9.8 

/six.  y.  :) 

=  -  r 

on  the  fin 

9.9 

■ r 

on  the  stabilizer 

9.10 

Figure  17  shows  the  panel  distribution  used  in  calculating  generalised  airforces  for  Staix’s 
T-tail.  In  order  to  highlight  the  panel  arrangement  in  the  vicinity  of  the  fin-stabiliser  junction 
only  the  starboard  half  of  the  stabiliser  has  been  included.  The  present  aerodynamic  ideali¬ 
sation  of  Stark's  T-tail  is  quite  similar  to  that  used  by  Kalman,  Rodden  and  Geising  [Ref.  12]. 
In  particular,  both  meshes  arc  refined  in  the  vicinity  of  the  fin-stabiliser  junction  and  the  tips  of 
the  stabiliser  in  order  to  improve  the  estimates  of  pressures  in  these  regions. 

In  developing  a  suitable  panel  arrangement  for  Stark's  T-tail  care  was  taken  to  ensure  that 
at  the  fin-stabiliser  junction  none  of  the  aerodynamic  panels  overlapped.  This  feature  of  the 
panel  distribution  is  clearly  evident  in  Figure  17  as  the  spanwise  lines  on  the  fin  and  stabiliser 
meet  at  the  fin-stabiliser  junction. 

If  the  panels  at  the  fin-stabiliser  junction  overlap  each  other  (i.e.  the  spanwise  lines  on  the 
fin  and  stabiliser  do  not  line  up  at  the  fin-stabiliser  junction)  it  has  been  found  that  the  pressure 
distributions  in  the  vicinity  of  the  junction  differ  from  those  obtained  for  the  lined  up  case. 
Figure  18  shows  plots  of  the  chordwise  variation  of  nondimensional  pressure  A  on  the  fin  and 
stabiliser  at  the  fin-stabiliser  junction  for  mode  /i  at  a  frequency  parameter  k  -  0  and  Mach 
number  .\t  0.  F-'or  the  case  where  the  panels  were  lined  up  at  the  junction  the  panel  distri¬ 

bution  of  Fig.  17  was  used.  Another  configuration  was  then  Investigated  where  the  panel  distri¬ 
bution  on  the  stabiliser  was  varied  from  that  of  Fig.  17  in  that  the  section  of  the  chord  from  the 
trailing  edge  of  the  fifth  chordwise  panel  to  the  trailing  edge  of  the  stabiliser  was  divided  into 
SIX  evenly  spaced  panels.  Only  the  rear  half  of  the  pressure  distribution  on  each  surface  is  shown 
since  the  effect  of  the  panel  overlap  was  localised  to  this  region. 

Figure  IX(a)  shows  Ifiat  the  pressures  on  the  fin  arc  only  slightly  modified  by  the  overlap 
of  panels  at  the  fin-stabiliser  junction.  In  contrast.  Figure  l8(/r)  shows  that  the  overlap  of  panels 
has  significantly  modified  the  pressure  distribution  acting  on  the  rear  half  of  the  stabiliser, 
although  only  two  pressure  points  have  been  significantly  changed. 

F  rom  the  above  example  it  is  evident  that  the  panels  should  be  made  to  line  up  at  any 
lunctions  if  accurate  anvl  cimsistent  predictions  of  tlie  pressure  distribution  are  desired,  f  or 
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this  particular  case  the  elTect  of  panel  overlap  did  not  modify  the  generalised  forces  to  any  signi- 
cant  degree  ;  the  variation  was  of  the  order  of  that  which  could  be  obtained  by  choosing  a  slightly 
different  (but  still  lined  up)  panel  distribution.  However,  for  an  alternative  set  of  modes  that  are 
more  sensitive  to  the  pressures  in  the  region  of  overlap,  the  elfect  could  have  been  significant, 

Table  5  presents  generalised  airforces  calculated  for  zero  frequency  (A-  0)  and  two  Mach 

numbers  M  0  and  M  0-8.  The  generalised  forces  ^n,  ().>i  and  (7ai  represent  the  yawing 
moment  due  to  yaw,  the  sideforce  due  to  yaw  and  the  rolling  moment  due  to  yaw.  The  agreement 
between  the  results  of  the  present  method  and  those  obtained  by  dilferent  workers  is  generally 
good.  If  we  focus  attention  on  the  doublet  lattice  calculations  of  Kalman,  Rodden  and  Geising 
[Ref  12],  it  is  seen  that  the  generalised  airforces  agree  well  with  the  largest  dilference  being  in 
Qn.  This  generalised  force  is  the  lin  yawing  moment  due  to  yaw.  and  since  the  side  force  acting 
on  the  fin,  ^ji,  agrees  well,  this  implies  that  while  the  present  method  and  that  of  Reference  12 
calculate  the  same  force  acting  on  the  fin  they  dilfer  in  the  chordwise  location  of  the  centre  of 
pressure. 

Table  6  presents  generalised  airforces  calculated  for  Mach  number  /V/  0-8  and  two 

frequency  parameters.  A  0  6  and  A  0'9.  The  results  obtained  by  Davies,  Zwaan  and 
Kalman.  Rodden  and  Cieising  were  obtained  from  Reference  12.  If  we  look  at  the  results  for 
A  0-6  and  M  0  8.  we  see  that  the  results  for  the  present  method  lie  well  within  the  range  of 
values  predicted  by  the  other  methods.  However,  the  agreement  between  the  magnitude  of  the 
generalised  forces  Du  and  Di-  calculated  by  Farrell  and  the  present  method  is  comparatively 
poor.  These  generalised  forces  represent  the  yawing  moment  due  to  yaw  and  sidesway,  and  are 
obtained  by  a  chordwise  weighting  I'f  the  pressure  distributions  due  to  modes  f]  and  fs.  Since 
Farrell  used  I'lily  four  panels  down  the  chord  of  the  fin  and  stabiliser  [Ref  2],  compared  to  nine 
and  eleven  for  the  present  method,  while  the  spanwise  panel  distributions  were  similar,  the  poorer 
definition  of  the  chordwise  pressure  distribution  thus  obtained  leads  to  quite  dilferent  estimates 
of  (All  and  (Au  c.ilculated  by  the  two  methods.  A  second  coarser  chordwise  panel  distribution 
was  used  with  the  present  method  with  the  result  thiil  the  generalised  forces  gn  and  were 
in  better  agreement  with  those  obtained  by  Farrell.  A  comparison  of  the  results  is  presented  in 
Table  7  where  it  is  seen  that  the  agreement  with  the  results  of  Farrell  is  improved  by  the  use  of 
the  coarser  chordwise  panei  distribution. 

Note  that  the  (  ase  2  panel  distribution  consisted  of  fitur  panels  down  the  chord  of  the  fin 
and  five  panels  down  the  chord  t'f  the  stabiliser.  This  combintition  was  chosen  to  allow  the  panels 
iilong  the  fm-siabihser  junction  to  line  up.  whereas  with  the  panel  distribution  used  by  Farrell 
this  would  nt't  have  been  possible,  fhis  feature  is  likely  to  h;ive  contributed  to  the  dilference 
between  the  Case  2  results  of  the  present  method  and  iho.se  of  Farrell. 


9..f  ONFR.V  Horizontal  Wing  and  I'ailplane 

file  (AM  RA  horizontal  wing  .iiul  tail  combination  is  shown  in  Figure  19.  In  this  particular 
evaniple  the  wing  and  tailplane  are  considered  to  be  identical  rectangles  of  chord  c  0  098 
metres  and  semi-span  v  0151.^  metres.  I  he  le.idmg  edge  of  the  tailplane  is  a  distance 
r.\  metres  downstream  of  the  trailing  edge  of  the  wing  aiul  the  plane  of  the  tailplane  is  a  distance 
(7/  metres  above  the  plane  of  the  wmg.  I  he  ptisition  of  the  tailplane  rehitive  to  that  of  the 
wing  IS  then  charticleriseil  by  the  noii'limenstonal  separation  partimeter  ,\  .iiul  height  partimeter 
II. 

In  Reference  \  Davies  h.is  presented  results  obttiincvl  using  his  lifting  surface  method  for  a 
set  of  fnir  combm.itions  of  freviuency  partimeter  and  Mach  number.  These  tire  given  below. 
V  iz. : 

\t  (1  fo  0  4:^  0  t.s  0  so 

A  I'  'ssp  It  24 II  I  s  I  .f  0-1112 

|S 


i 

The  frequency  parameter  k  is  based  on  a  reference  length  equal  to  the  chord  of  the  wing.  The 
variation  in  Mach  number  and  frequency  parameter  corresponds  to  an  oscillation  of  almost 
constant  angular  frequency  in  a  compressible  flow  of  varying  velocity. 

The  modes  of  oscillation  that  were  studied  consisted  of  four  rigid  body  modes,  namely 
i  vertical  heave  of  the  wing,  vertical  heave  of  the  tailplanc,  nose-up  pitch  of  the  wing  about  its 

I  mid-chord  and  nose-up  pitch  of  the  tailplanc  about  its  mid-chord.  These  four  modes  can  be 

described  analytically  as 


./i(-v,  r,  --) 

1  ■() 

on  the  wing 

9.11 

0 

on  the  tailplanc 

9.12 

,1'Ax.  .f.  --) 

0 

on  the  wing 

9.13 

1  ■() 

on  the  tailplanc 

9.14 

y,  :) 

.V  1 

<  2 

on  the  wing 

9.15 

0 

on  the  tailplanc 

9.16 

l.{.\.  A'.  --) 

1' 

on  the  wing 

9.17 

on  the  tailplanc 

9.18 

w  here  the  origin  of  the  coordinates  is  at  the  centre  of  the  leading  edge  of  the  w  ing. 

In  calcukitit'ns  with  the  present  method  a  panel  distribution  of  ten  evenly  spaced  chord- 
wise  panels  ;ind  twelve  evenly  spticed  spatiwise  panels  was  used  for  both  the  wing  and  tailplanc. 

f  rom  the  evteiisive  range  of  results  obtained  by  Dtivies.  which  cover  a  range  of  values  of 
.\  from  0  to  A  -/  for  //  0  and  II  IS.  the  nonplanar  case  [H  IS)  with  A  2 

and  I  2  were  chosen  to  be  used  in  ctilculations  with  the  present  metluvd.  For  this  choice  the 
wing  tailplanc  provimity  elfects  wvvuld  be  dilferent  from  those  encountered  for  the  very  closely 
spticed  .ACi/XKD  wing  and  ttiilplane  configur.ttion,  thus  prvniding  ;i  somewhat  dilferent  situation 
compared  to  the  ,\(i,-\Rl)  conftguration. 

Tiible  S  gives  results  for  the  v.irious  \l  k  comhintitions  ;is  obtained  by  Dav  les  and  the  present 
method  for //  I  S  and  A  2.  C  vimparison  of  the  twiv  methods  indicates  th;it  there  is  generally 
good  agreement  in  the  m.ignitiKles  of  the  genertilised  forces,  the  largest  relative  dilferences 
occurring  in  the  terms  0\i  and  (J-.n.  These  represent  induced  forces  and  pitching  moments 
on  the  wing  due  tv>  oscilhitii'n  id  the  ttiilpl.ine.  ;ind  ;ls  such  they  are  more  likely  to  vary  since  for 
this  pttrlicuhir  case  they  represent  a  sccoiul  v'rder  etfect. 

Turning  our  attention  to  the  c.ilctilated  phtise  .ingles  it  is  seen  that  there  is  good  agreement 
at  the  highest  Mach  nutnher.  1/  O  N.  I  hen.  as  the  Mach  number  decreases  the  tigreement  gets 
progressively  worse  until  ,it  the  lowc't  Mach  number  the  ililference  in  ph.tse  tingles  falls  between 
3  tind  5  degrees.  Ibwvevci.  this  is  only  mtirgintilly  wnrse  tluin  the  mtivimum  dilTerence  of  I  -5 
degrees  encounterevl  for  the  highest  Vl.ich  number. 

l.ible  d  gives  results  .IS  obt.iineil  by  I). ivies  .nul  the  present  method  for //  I  ,N  and  I  2. 
It  IS  seen  th.it  the  agreement  m  the  genertiliscil  forces  (>|i  and  (fti  htis  been  improved  comptircd 
with  the  results  shown  in  l.ible  S  It  is  cle.ir  th.it  by  reducing  the  induced  forces  and  pitching 
motnents  on  the  wing  iluc  to  oscill.ition  of  the  tailplanc  have  increased. 


in.  INTI  RAt  tivt;  (;k\i>mk  s  disim  ax  i  aiim  ia 


When  devell'p'ng  ,i  d.ii.i  lile  fur  the  ulc.ihstitinn  of  .i  given  ct'iiligiiration  it  is  desirable  to 
have  .ivail.ible  the  v.t|i,ibilily  to  genci.ite  ,i  plot  of  the  p.inel  distribution  corresponding  to  the 


input  data.  This  facilitates  the  checking  of  the  panel  distribution  for  an>  errors,  as  \sell  as  per¬ 
mitting  the  visual  comparison  of  various  alternatives. 

The  present  program  incorporates  the  ability  to  generate  a  data  file  from  a  goen  panel 
distribution  which  is  compatible  with  the  PAFEC  Unite  element  package  that  is  m  common 
use  at  ARL.  This  file  is  then  processed  by  PAFEC  so  that  the  PAl  Ft  Interactoc  'ir  phics 
System  (PIGS)  can  be  used  to  inspect  the  panel  distribution.  PKiS  pnnides  mans  powerful 
commands  that  allow  the  user  to  study  the  panel  distribution  in  great  detail  from  aiu  chosen 
angle.  A  window  facility  permits  close  inspection  of  junctions  at  the  intersections  of  surfaces, 
thus  these  critical  areas  can  be  easily  checked  to  ensure  that  panels  line  up  correctK 

Very  complex  panel  distributions,  such  as  those  in\ol\ing  multiple  underwmg  stores  with 
aerodynamic  surfaces,  can  be  checked  with  ease.  This,  together  with  the  abilitc  to  obtain  hard 
copies  of  any  plots  generated  at  the  VDU  terminal,  makes  PICJS  an  indispens.ihle  tool  for  the 
validation  of  panel  distributions. 

The  inherent  compatibility  of  a  doublet  lattice  panel  distribution  with  mi'st  finite  element 
graphics  packages  means  that  the  job  of  preparing  a  data  file  for  such  a  package  is  a  relatively 
simple  task.  In  program  AIRFORTP  only  one  additional  subroutine  was  required  and  only 
this  subroutine  would  need  to  be  modified  if  a  finite  element  graphics  package  other  than  PIGS 
were  to  be  used.  The  advantage  of  such  an  approach  is  that  existing  computer  graphics  systems 
can  be  utilised,  thus  obv  iating  the  need  to  develop  additional  programs  to  improve  and  speed  up 
the  data  development  and  checking  phases  involved  in  defining  a  panel  distribution. 

The  possibility  of  creating  the  entire  panel  distribution  using  PIGS  was  considered,  only 
to  be  rejected  for  a  number  of  reasons.  Firstly,  the  common  requirement  that  panels  on  different 
surfaces  line  up  streamwise  impo.ses  significant  constraints  on  the  panel  generation  scheme  which 
could  lead  to  difficulties.  Secondly,  certain  types  of  panel  distributions  cannot  he  developed 
easily  using  PIGS,  since  there  may  be  many  geometrical  considerations  that  must  be  taken  into 
account.  The  present  method  which  utilises  a  data  input  technique  tailored  to  the  requirements 
of  program  AIRFORTP,  as  well  as  employing  the  graphics  capabilities  of  PIGS,  provides 
capable  and  easy  to  use  data  preparation  techniques. 


11.  DISCUSSION 

The  doublet  lattice  method  is  a  simple,  versatile  and  accurate  liffing  surface  theory  which 
is  capable  of  analysing  aerodynamic  surfaces  of  arbitrary  planform  and  dihedral.  Although 
they  were  not  considered  in  this  report,  control  surfaces,  either  full-  or  partial-span  may  be 
included.  Problems  of  intersecting  and  interfering  nonplanar  configurations,  such  as  a  vving- 
pylon-store  combination,  a  T-tail  or  V-tail.  a  wing-tail  combination,  etc.,  may  be  analysed. 

The  computer  program  AIRFORTP  can  handle  most  configurations  provided  that  they  are 
idealised  within  the  restrictions  imposed  by  the  doublet  lattice  method.  Hence  the  previous 
capability  at  ARL  for  analysing  vving-tail-lin  combinations  has  been  significantly  extended. 
Appendix  Ci  contains  general  considerations  in  regard  to  the  development  of  panel  distributions 
used  to  model  any  given  aerodynamic  configurations.  Appendix  H  gives  a  description  of  the 
format  of  the  input  data  file  and  mode  generation  subroutine  used  by  program  AIRFORTP. 
An  example  of  an  input  data  file  and  mode  generation  subroutine  is  given  m  Appendix  I  for  the 
analysis  involving  Stark  s  F-tail  configuration. 

The  results  obtained  by  the  preseni  method  for  planar,  nonpl.inar  and  intersecting  lifting 
surface  combinations  have  been  compared  with  the  results  of  dilVerent  workers,  fhe  latter  were 
obtained  by  both  doublet  lattice  ;md  lifting  surlace  methods  fhe  results  of  the  presetit  method 
are  seen  to  lie  within  the  range  of  results  obtained  by  the  other  methods.  Problems  with  a  sintill 
vertictil  (non-/ero)  septiration  between  streamwise  columns  of  panels  ean  be  handled  with  ease 
up  to  the  point  where  the  non-zero  separation  is  so  small  as  to  he  negligible. 
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12.  CONCLUSIONS 


The  doublet  lattice  method  is  readily  applicable  to  a  large  class  of  interfering  lifting  surfaces. 
As  currently  programmed  at  ARL  the  method  can  be  used  for  general  nonplanar  and  nonparallel 
lifting  surfaces,  both  intersecting  and  nonintersecling  in  nature.  Results  have  been  obtained  for  a 
number  of  configurations  and  comparisons  of  the  results  have  been  made  with  results  obtained 
by  other  workers.  The  comparisons  show  generally  good  agreement. 
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APPENDIX  A 


Definitions  relating  to  a  general  combination  of  a  sending  and  receiving  panel 

Let  us  consider  the  interaction  of  a  receiving  panel  r  and  a  sending  panel  s  which  are  defined 
in  a  global  axes  system  (A",  Y,  Z)  as  shown  in  Figure  2.  These  panels  are  inclined  with  dihedral 
angles  defined  as  yr  and  ys-  The  coordinates  of  points  A  and  B  at  either  end  of  the  doublet 
line  on  the  sending  panel  are  A  :  (xa,  >’a,  ^a)  and  B  :  (xb,  >'b,  2b).  The  coordinates  of  the  control 
point  C  on  the  receiving  panel  are  C ;  (x,  y,  z). 

The  position  vector  corresponding  to  points  along  the  doublet  line  may  be  written  in  terms 
of  the  variable  t  as  follows : 


rAn(t)—  .VAB(/)i+  Tab!/)] +  ZAn(t )k  0  <  /  ^  1  A.l  | 

( 

where  i 

A'ab(0  =  JfA  +  (.rB-XA)  /  A. 2 

Tab(/)  =  .i'A  +  (  Vb— J’a)  t  A. 3  I 

Z\n(i)  =  za  +  (zb”Za)/  a. 4 


From  the  definitions  of  X|,  yi  and  zi  in  Equation  4.1  it  is  clear  that  when  we  calculate  the 
line  integral  of  Equation  3.2  the  variables  .vi,  i  i  and  zi  may  be  written  in  the  following  form  : 

xi  =  a^  +  bxi 

>•1  =  Oy  +  byt  A. 5 

2i  =  az  +  bit 

where 


Ux  =  X-Xa 

fly 

=  r-jA 

flz  =  2  — Za 

A.6 

Ax  .Va-  Xb 

by 

-  Va  ”3’1! 

1! 

N 

> 

1 

A. 7 

In  the  equations  for  the  Kernel  function  the  variable  r'^  is  present.  The  equation  for  r’  may 
be  written  as : 


where 


r- 

.n-  +  2i^ 

flrt  t  brt  <'r 

A.8 

Ur 

Ay-+6,- 

A.9 

hr 

2[fli/)v  i  fl/A/) 

A. 10 

(  r 

fly-  ->  It,  ' 

A.l  1 

If  the  minimum  value  of  n  for  a  given  panel  combination  is  denoted  by  and  the  location 
of  this  minimum  by  /min,  then  we  have  that ; 


train  —  —bxjlxit 


A.12 


r^rain  —  Cr  —  -- 

4ar 


A.13 


A  decision  which  determines  the  integration  procedure  to  be  used  is  made  on  the  basis  of  the 
values  of  /•-rain  and  /min  obtained  for  a  given  panel  combination.  Due  to  the  organisation  of 
panels  into  streamwise  columns,  /mm  and  r^min  need  only  be  calculated  for  combinations  of 
columns  rather  than  individual  combinations  of  panels. 


APPENDIX  B 


Calculation  of  the  infinite  integrals  Ii  and  I2 

In  calculating  the  Kernel  function  two  infinite  integrals,  h  and  /2,  occur.  For  completeness 
the  expressions  are  repeated  below : 


h{ui,  v) 


U+M-)®’ 


du 


B.l 


/2(M1,  1/) 


00 


Ml 


(1+1/2)5  2 


.du 


B.2 


where  the  symbols  mi  and  v  were  defined  in  Section  4.1. 

From  Reference  7,  integration  of  h  by  parts  once  and  integration  of  3/2  by  parts  twice 
leads  to : 


where 


3/2(Mi,  .)  =  C->--.|(2  +  /Vm.)[i  2]-(,  +r,2)3  2 

— /V/o(Mi,  i’)  +  v'-Jo(Ul,  v)! 


,/,  L' 


X 

u 

c''  "1 

u 

(1  +m2)I  2 

Hi 

B.3 


B.4 


B.5 


B.6 


The  integral?.  A.  and  ./n  can  be  e\aluated  using  approximations  to  m(I+m2)  ’  -  developed 
b\  Laschka  [Ref.  S|  in  an  exponential  form  for  11  ^0: 


II 

(I  +m2)> 


II 


u  I 


B.7 


where  r  0  -372  and  the  //„  are  given  in  Table  10. 


The  integrals  /o  and  /o  for  i/i  >  0  then  become ; 


/o(M1,  •') 


11 


(n^c^  +  v^) 


(nc—iv) 


n  =  1 


B.8 


•/o(Wl,  v) 


11 


cine  , 

,  2  2.  '2\2^''*‘’^ ~v^  +  ncuiin^c^  + 1'^) 


n  =  1 


B.9 


~iv[2nc  +  ui(n'^c'^  +  v^)]} 

Symmetry  properties  of  the  integrands  of  h  and  I2  have  permitted  the  consideration  of 
only  non-negative  arguments  ui  since,  for  wi  <  0,  we  have  that : 

/i(wi.  v)  =  2/Je{/i(0, .)}  -/?e{/i(-i/i,  v)}-(-;74/i(-wi,  v)}  B.IO 

h{uu  v)  =  2/?e{/2(0,  k)} - Re{I-^ -u,  i^)} -f- i/miH-u,  v)}  B.  1 1 

Hence  the  desired  integrals  /i  and  h  have  been  explicitly  defined  for  the  range  of  cases 
ui  ^  0  and  wi  <  0. 


APPENDIX  C 


Integration  procedure  when  r^  >  0  over  the  integration  interval 

The  normal  wash  influence  coefficient  relating  the  normal  wash  on  the  Tth  panel  to  the  pressure 
over  the  j'th  panel  is: 


Dii 


Cj  COS 
47r 


KiidI 


h 


C.l 


where  /--  <  0  over  the  interval  /j.  By  analytically  incorporating  the  \  r~  behaviour  into  the  in¬ 
tegral,  Equation  C.  1  becomes : 


flij  - 


fj  cos  i3j 

47r 


di 


C.2 


where 


6(.vi,  >1.  n)  =  (KiTi  + 


C.3 


By  subdividing  the  interval  Ij  into  n  smaller  intervals  we  may  write: 


Dii  - 


cj  cos 
4jt 


^i  V 


m  =  I 


dl 


C.4 


where  the  line  integral  is  evaluated  over  w  straight  line  segments  between  points  Am:  (VAm, 
,i'Am.  CAni)  and  B,n :  (.vnm,  .Him,  Cum).  For  convenience  the  m.  i  and  ./  subscripts  will  be  dispensed 
with,  and  so  the  integral  of  interest  in  Equation  C.4  may  be  written  as : 


/ 


r- 


vj 


C.5 


Let  the  line  between  points  A  ;  (  va.  i  a.  z,\)  and  B  :  ( Vn.  i  n.  rn)  be  represented  by  the  para¬ 
metric  relationship: 


rAiil/’)  A  ah(/’)' +  )'An(/>)j  +  ZAU(/')k 


1  ^  ^  +  I  C.6 


where 


A  ah(/7)  V  f  \p  V  (aa  +  A  id  2 


V  (An  a-a)/2  C.l 


Ykb,{p)  =  y+yp  y  =  0'a-I->'b)/2 

y  =  0'b->’a)/2 

C.8 

Zab(/>)  =  z  +  zp  f  =  (za  -1-  zb)/2 

n  =  (zb-Za)/2 

C.9 

Equation  C.5  may  now  be  written  in  the  form  : 

^  '  r G(xi{p),  yi(p),  zi(p))  dl 

_i  3’i(p)*  +  zi(p)“  dp  ^ 

*  % 

CIO 

where 

xi(p)  =  .vo-A'ab(p) 

C.ll 

yi{p)  =  yo—  yAB(p) 

C.12 

2\(p)  =  Zo-Zab(/>) 

C.13 

and  (.Vo,  >’o,  zo)  is  the  location  of  the  control  point  for  the  panel  combination  being  considered. 
Now  we  have  that : 

dl  .  .  , 

y  —  (Tab  •  fAfl)'  ■ 

dp 

=  (Xab(p)^+  Yab(p)^+Zab(pW'^ 

which  yields  the  result : 

dl 

—  =  {x^+y^  +  z^y^ 
dp 

C.14 

By  substituting  Equations  C.ll,  C.12  and  C.I3  into  Equation  C.IO  we  obtain: 

dp_,  D(p)  ^ 

*  % 

C.15 

where 

D(p)  =  ap^  +  bp+c 

C.16 

a  —  y^  +  P 

C.I7 

b=  -2[)’{yo-y)  +  z(zo-z}] 

C.18 

c  =  (3'o~/)2-I-(zo-zF 

C.19 

Now,  in  order  to  evaluate  the  integral  in  Equation  C.I5,  let  us  represent  G{p)  by  a  second-order 
Lagrangian  interpolation  polynomial : 


3 

(aip^+b\p+  ci)G(p\) 


C.20 


where  the  abscissas  pi  arc  the  zeros  of  the  Legendre  polynomial.  Substituting  Equation  C.20 
into  Equation  C.15  yields  : 


I 


3 


C.21 


where 


W, 


+  1 

a, 

I 


P- 

D 


(Ip  +  b, 


4/ 


+  1 


t'i 


C.22 


From  the  tables  of  integrals  presented  b>  Dwighi  in  Reference  14  we  lind  the  following  indefinite 
integrals  corresponding  to  the  definite  integrals  in  Equation  C',22  : 


p  h  h-  2(ic 

(I  2a-  2a- 


D 


dp 


C.23 


P  .  \  b  I 

dp  --  In  />  -  dp 

n  2a  n  ' 


C.24 


I  .  2  b+2ap 

-  dp  tan  ' 

D  (4ac  -  b-)' -  (4«(  -6-)' 


.  b‘^~4ac<0  C.25 


I  I  ^  ^  ■ 

(b-  4at}‘-  b-^2apr(b-  4(;( )' - 


.  b-  4ac  >  0 


2 

h  +  2ap 


,  b-  -  4ai  0 


The  coelhcionts  (/,.  and  ci  of  Equation  C  .20  may  be  calculated  from  : 


a\p-(h\p^c\  (p  p>){p  />:))(/)!  p-){p]  /13)  C.26 

a2p-  +  h:p^c>  ip  pi)(p  p:t)  (p2  p\){pi  p-.d  C.27 

a-.\p-  *  h^p  (p  p\){p  pd  ip.\  P\){p-.i  Pi)  C.23 

b\  equating  like  powers  of />.  Thus  we  obtain  : 

</i  I  (/>!  /’jt(/»i  P:\)  h\  (/>_•  t /»:i)ui  (  |  Pip^cp  C.29 

a-i  I  (/’_•  /ii)(/>-  p:d  b:  {p\>pdai  i  j  p\PMii  C..30 

i/:i  I  (/t:i  /’iM/»3  Pi)  b-.t  (Pi‘Pi)a:\  ( :i  />i />_■</;)  C.3I 


where  the  abscissas /»,  and  coetlicients  a,,  hi  and  <  ^  are  presenteil  in  fable  I  I. 


I 


By  choosing  that  the  abscissas  pi  correspond  to  the  abscissas  in  the  common  Gaussian 
quadrature  formula : 


Ap)dp 


it  is  anticipated  that  the  integration  accuracy  will  be  maximised  compared  to  any  other  choice 
of  Pi. 

Note  that  since  is  greater  than  or  equal  to  zero  then  the  determinant  b'^—^ac  is  less  than 
or  equal  to  zero.  Hence  in  Equation  C.25  the  integral  when  b'^  Aac  >  0  is  not  applicable  to 
the  problem  at  hand  and  should  be  neglected. 


APPENDIX  D 


Integration  procedure  when  r-  0  in  the  centre  of  the  integration  inter\al 

in  the  doublet  lattice  method,  when  the  r-  0  smgulantv  i><.curs  it  must  do  so  in  the  centre 
of  the  interval  /j  over  which  the  line  integral  of  Equation  is  dclined.  f  ollowmg  the  method  of 
Appendix  C,  the  interval  I,  may  be  subdivided  into  ii  smaller  intervals.  Eurthermore,  if  we  assume 
that  the  singularity  falls  symmetrically  vvithin  the  A'th  interval,  then  the  normalwash  influence 
coefficient  Dij  may  be  written  as  ; 


Cl  cos  dj( 
477  I 


A  I 

X 


D.l 


n 

\ 

*T+ 1 


B.. 

A„ 


The  integrals  in  the  first  and  third  terms  may  be  calculated  using  the  method  detailed  in  Appendix 
C.  It  now  remains  to  derive  a  method  for  solving  the  finite  part  integral  in  Equation  D.l.  If  we 
dispense  with  the  subscripts  then  the  integral  of  interest  is  : 


/  - 


D.2 


Now  since  r-  0  in  the  centre  of  the  interval  from  point  A  ;  (.v a.  .I'a.  to  point  B  :  (  vn. 
I  n,  rn).  the  line  between  these  two  points  can  be  represented  by  the  parametric  relationship : 


r(/7) 

.v(/))i  +  i  (/i)j4  :(r)k 

1 

$  ^  4  1 

D..^ 

Mr) 

.V  ^  V/7 

V  (.Ylt 

•Va)  2 

V  (.Va4V| 

;l  2 

1)4 

Mr) 

In  +  177 

4  (I'll 

I'a)  2 

Dr 

--(r) 

rn  4  f/> 

f  (rn 

r.A)  2 

17.6 

and  the  control  point  is  located  at  (  Vn,  in.  rn).  Note  that  when  r-  0  the  eoiiiri'l  point  .iiui  lift 
point  (shown  in  Eia.  2)  lie  on  the  same  streamwise  line. 


Using  the  above  definitions  of  .v(/>),  >’(p)  and  z(p)  we  obtain  ■. 
xi(p)  =  Aro-(.v  +  xp) 


D.7 


yi(p)  =  -yp 

z\(p)  =  -  zp 

Hence  we  see  that  r’  =  yi^  +  z\'^  is  simply  : 

,2  (yi+:Z)pi 


and,  as  in  Appendix  C,  we  have  that : 


=-  (.xi+F+«)'  ^ 


Thus,  using  Equations  D.7  to  D.ll,  Equation  D.2  becomes: 


'  r  K‘* 


Following  the  method  of  Farrell  [Ref.  2],  the  numerator  G(p)  may  be  approximated  by  a 
polynomial  of  order  (n  I ) . 


Q(p)  ^  hypi  'C(/?i)  D.I3 

where  the  hy  represent  the  coefficients  of  the  polynomial.  The  abscissas  Pi  are  the  roots  of  the 
/;'th  order  Tchebycheff  polynomial  of  the  first  kind.  The  abscissas  may  be  determined  from  the 
following  equation  [Ref.  15.  p.  889): 


p,  cos 


( 2/  1  )7r 


By  substituting  Equation  D.l.''  into  Equation  D.I2  we  obtain  : 


I  <//  \‘ 

/  ^  r,6(/8) 

r-  +  c-  (Ip  , 

/  I 


.  4  I  /• 


Xpj  '-'(ip. 


By  equating  like  powers  of  p,  the  Aji  may  be  calculated  from : 


V 


where  11  is  the  restricted  product  for  «  =  1,  2 . « ;  u  ^  i, 

u 

and  n  is  the  restricted  product  for  r  =  1,2 . n  ;  r  #  /. 

V 


As  in  Reference  2,  an  eighth  order  polynomial  has  been  chosen.  This  corresponds  to  n 
and  the  abscissas  pi  and  weights  Ci  have  been  calculated  and  are  listed  in  Table  12. 


APPENDIX  E 


Estimating  the  position  of  the  zero  crossings  of  Kor^ 

Let  us  consider  the  interaction  of  a  receiving  panel  r  and  a  sending  panel  s.  We  wish  to 
determine  estimates  of  the  positions  of  the  zero  crossings  of  the  numerator  of  Equation  5. 1  for 
the  steady  case.  These  zero  crossings  will  occur  somewhere  along  the  line  which  is  colinear 
with  the  doublet  line  in  the  sending  panel,  and  we  are  specifically  interested  in  the  case  where 
they  occur  on  the  doublet  line  itself. 

The  numerator  discussed  above  is  defined  as : 

Aor‘^  =  cos  (yr— ■>'s)+  cos  yr— .Vi  sin  yrM^i  cos  ys  — >'i  sin  ys)lf^  E.  1 

Multiplying  through  by 

=;.,2  +  Zi'^ 

we  obtain  : 

AoH  cos  (yr  ysK.i'i-  +  ri-)  +  cos  yr  -  >1  sin  yr)(ri  cos  ys  -  Ji  sin  ys) 


E.2 

By  substituting  for  n  and  ,i  i  using  Equations  A. 5  to  A.l  1  (see  Appendix  A)  it  is  possible  to  write 
Equation  E.2  in  the  following  form  : 


Kor^  =  Ar^  +  Bt-\-C  E.3 

where 

A  =  A-i'«'[cos  (yr  -ys)flr)  +  A:2's>[a2a4]  E.4 

B  =  Ai''’’[cos(yr-  ys)i>r]  +  /k2'^’[ata4  +  n2<J3]  E.5 

C  -■=  Ai'-^'fcos (yr- -  ys)cr]  +  /Ci's'ia.fls]  E.6 

ar  =  6y2  +  ^z-  E.7 

bi  =  2[ayby  +  a,.b,]  E.8 

cr  =  ay-  +  a/  E.9 

fli  --f/zcosyr  r/ysinyr  E.IO 

«2  /’zCOSyr  sin  yr  E.ll 

«3  «z  cos  y^  f/y  sin  y.  E.  1 2 

(i\  by  cos  /7y  sin  y.  E.1.3 


In  the  above  equations  it  is  assumed  that  and  are  constant  for  the  panel  com¬ 
bination  under  consideration. 


Inspection  of  Equation  E.19  leads  to  the  observation  that  there  is  only  one  zero  crossing 
when  the  sending  and  receiving  panels  are  perpendicular  to  each  other.  An  example  of  such 
behaviour  is  presented  in  Figure  5. 

The  positions  of  the  zero  crossings  (in  terms  of  the  parameter  l)  of  the  numerator  of  Equa¬ 
tion  5.1  can  be  estimated  by  solving  the  following  equation  : 

At2  +  Bt  +  C  =  0  E.26 

for  the  values  of  I. 

In  the  above  analysis  it  has  been  assumed  that  and  may  be  taken  to  be  constants 
over  the  length  of  the  doublet  line  in  the  sending  panel.  This  approximation  is  valid  since 
both  K\'^'  and  are  slowly  varying  functions  and  the  interval  over  which  they  are  assumed 
to  be  constant  is  itself  quite  small.  More  often  than  not,  the  zero  crossings  will  occur  outside 
the  interval  defined  by  the  doublet  line.  In  this  case  the  assumed  values  of  A'l'*’’  and  A:'2“’’  may 
not  be  accurate,  but  this  is  unimportant  since  we  are  only  interested  when  the  zero  crossings 
occur  along  the  doublet  line.  When  the  zero  crossings  occur  along  the  doublet  line  the  estimate 
of  their  location  may  be  improved  by  repeating  the  calculation  using  revised  average  values  of 
and  A'2'*’  calculated  over  a  smaller  interval  defined  by  the  initial  estimates  of  the  zero 
crossings. 

For  critical  cases  where  the  zero  crossings  are  close  together  the  approximation  that 
and  are  constant  over  the  small  interval  is  particularly  good.  Geising  et  al.  [Ref.  7]  have 
shown  that  when  /■--‘O  the  terms  Ki-*2  and  where  Ki  and  K2  are  elements  of  the  un¬ 
steady  Kernel  function  given  in  Equation  4.1.  Hence,  as  r^-*0  then  and  4. 


APPENDIX  F 


Estimating  the  position  of  the  turning  points  of  Kor-  for  a  perpendicular  combination  of  sending  and 

receiving  panels 


Let  us  consider  the  interaction  of  a  perpendicular  combination  of  receiving  panel  r  and 
sending  panel  5.  Utilising  the  expressions  developed  in  Appendix  E,  the  numerator  of  Equation 
5. 1  can  be  written,  for  the  steady  case,  as : 


= 


Afi  +  Bt  +  C 

Orf^  +  brf  +  Ct 


F.l 


where  the  coefficients  A,  B  and  C  are  defined  in  Equations  E.19  to  E.2I  and  the  coefficients  r/r. 
Arand  cr  are  defined  in  Equations  E.7  to  E.9.  Since  the  panels  are  perpendicular  cos  (yr  ~  ys)  -  0, 
and  Equation  F.l  becomes; 


+  a\) 

Aor-  =  '  r.Z 

4-  b[t  4“  Ct 

where  the  coefficients  ai,  a-i  and  as  have  been  defined  in  Equations  E.22  to  E.24. 

If  we  assume  that  is  constant  for  the  panel  combination  under  consideration,  then  : 

d(K(\r'’)  a3K>*^'{ata>t'-  +  2arair  +  btai—CTai) 

dt  {ati“  '^bft -c  Ct)'^ 

After  setting  the  above  expression  to  zero,  the  positions  of  the  turning  points  may  be  estimated 
by  solving  the  following  equation  : 


ara>l  -  +  2arai  /  +  brai  —  (ra>  =  0 


F.4 


for  the  values  of  t. 


APPENDIX  G 


Considerations  in  the  selection  of  panel  distributions 

In  the  doublet  lattice  method  the  surface  is  divided  into  a  number  of  panels.  A  number  of 
rules  and  considerations  apply  to  the  choice  of  panel  distribution  for  a  given  configuration ; 
these  are  as  follows : 

1.  The  panels  are  trapezial  in  shape  and  the  sides  of  panels  must  be  aligned  streamwise. 
Panels  are  arranged  in  strips  parallel  to  the  free-stream.  The  aspect  ratio  of  the  panels 
should  not  be  large.  For  the  unsteady  case,  an  aspect  ratio  of  order  unity  or  less  is 
preferred  [Ref.  7]. 

2.  Surface  intersections,  surface  edges,  control  surface  boundaries,  fold  lines  and  other 
discontinuities  should  lie  on  panel  boundaries. 

3.  For  wing-tail  type  configurations  the  strip  boundaries  on  the  tail  must  be  aligned  with 
those  on  the  wing. 

4.  The  analysis  assumes  that  the  pressure  is  constant  over  each  panel.  Thus  the  results 
will  be  more  accurate  if  the  panels  are  smaller  where  the  pressure  distribution  varies 
more  rapidly.  Consequently  the  panels  should  be  concentrated  near  the  wing  tips,  the 
leading  and  trailing  edges  of  the  surface  and  near  control  surface  edges.  This  also 
applies  to  all  regions  where  the  normalwash  boundary  conditions  are  diseontinuous. 

5.  Where  surfaces  intersect,  the  panels  should  be  arranged  such  that  the  sides  of  panels 
situated  along  the  line  of  intersection  do  not  overlap.  This  means  that  the  spanwise 
lines  must  coincide  as  shown  in  Figure  17. 


APPENDIX  H 


Input  data  for  program  AIRFORTP 

Program  AIRFORTP  calculates  the  pressure  distribution  and  generalised  airforces  for 
general  lifting  surface  configurations  oscillating  in  subsonic  flow.  The  data  file,  AIRFOR  .  DAT, 
for  this  program  has  the  following  form  ; 

GNU,  AM,  RL 
IPLOT,  IPRNTP 
NONODE 

I,  XN(I),  YN(I),  ZN(1) 


NGRPS 

GROUPTITLE 

J,  JINODE(J),  J2NODE(J),J3NODE(J).  J4NODE(J) 
IREF(J,  1),  IREF  (J,  2),  IREF  (J,  3).  IREF  (J.  4) 
NC(J),  NS(J),  NCEVEN,  NSEVEN,  ITCORD 
(TCORDl  (J,  K),  K  -  1,  NC(J)+  I) 

(TCORD2  (J,  K),  K  1,  NC(J)+  1) 

(TSORDS  (J.  K),  K  l,NS(J)+l) 


NGRPSA 

(IGRPSA  (K),  K  -  1.  NGRPSA) 
NMODEA 

(IMODEA  (K),  K  ^  1.  NMODEA) 


where 


wl 

GNU  frequency  parameter, 

am  Mach  number.  AM  <10. 

RL  reference  length,  /. 

IPLOT  Flag  for  initiating  preparation  of  data  file  for  use  in  PIGS  graphics 

package. 

IPLOT  0  :  no  data  preparation  ; 

IPLOT  -  I  :  PAFEC  data  file,  BOX  .  DAT,  is  created  and 
program  AIRFORTP  stops  execution. 


IPRNT 


Flag  for  choosing  format  of  print  out  of  the  calculated  pressure 
distributions. 


IPRNT  =  0 
IPRNT  =  1 
IPRNT  -  2 


NONODE 


1 

XN(1) 

YN(1) 

ZN(I) 

NGRPS 


GROUPTITLE 


J 


JINODE(J) 

J2NODE(J) 

J3NODE(J) 

J4NODE(J) 

IREF  (J,  I) 
IREF(J,  2) 
IREF  (J,  3) 
IREF  (J.  4) 

NC(J) 

NS(J) 

NCEVEN 


NSEVEN 


ITCORD 


no  pressures  printed ; 

pressures  for  each  mode  printed  in  column  fashion  ; 
pressures  for  each  mode  printed  in  matrix  array 
form  for  each  group  of  panels. 

number  of  nodes  to  be  defined.  Nodes  need  not  be  numbered 
sequentially.  I  ^  NONODE  <  100. 

node  identification  number.  1  ^  I  <  100. 

X  coordinate  of  node  I. 

Y  coordinate  of  node  1. 

Z  coordinate  of  node  I. 

number  of  groups  (of  panels)  to  be  defined.  Groups  need  not  be 
numbered  sequentially.  1  $  NGRPS  ^  50. 

title  to  describe  the  group  whose  input  data  follows.  Up  to  70 
characters  are  allowed,  e.g.  PORT  WINGLET. 

group  identification  number.  1  ^  J  NGRPS.  Groups  need  not 
be  numbered  sequentially. 

Node  numbers  defining  the  leading  and  trailing  edges  of  group  J. 
The  sequence  I.  2,  3,  4  follows  the  group  node  numbering  con¬ 
vention  (defined  in  Section  8.1  and  Fig.  8). 


numbers  defining  the  senses  of  positive  pressures  on  the  source 
group  and  any  images  that  have  been  created.  They  may  take  the 
values  0.  i.  +  I  except  for  the  source  group  IREF  (J,  1).  which 
must  be  non-/.ero.  Sec  Section  8.3  for  a  detailed  description. 

number  of  chordwise  panels  on  group  J. 
number  of  spanwise  panels  on  group  J. 

flag  for  choosing  whether  or  not  panels  down  the  chord  of  group  J 
are  to  be  evenly  spaced. 

NC  EVEN  0  :  panels  not  evenly  spaced. 

NC'EVEN  I  :  panels  arc  evenly  spaced. 

flag  for  choosing  whether  or  not  spanwise  columns  of  panels  on 
group  J  are  to  be  evenly  spaced. 

NSFiVEN  0:  columns  nv)t  evenly  spaced. 

NSEVEN  I  :  columns  are  evenly  spaced, 

flag  for  choosing  whether  or  not  the  ordinates  of  chordwise  edges 
of  panels  for  the  side  containing  nodes  J2NODE(J)  and  J4- 
NGDE(.I).  rCORD2(J.K ).  are  to  be  provided  as  input  data. 
This  feature  allows  panels  of  constant  chord  to  be  located  on  a 
tapercil  lifting  surface.  If  ITCORD  0  then  the  line  of  data 

containing  the  ordinates  TCORD2(J.K)  is  omitted  and  the  pro¬ 

gram  defaults  to  TC()RD2(J.K )  TCORDKJ.K).  If  IT(  ORD 
I  then  the  line  of  data  containing  TCORD2(J.K)  must  be  in¬ 
cluded. 


hal^ 


1 


TC0RD2(J,K) 


TSORDS(J,K) 


NGRPSA 


TCORDl(J,K)  ordinates  of  the  chordwise  edges  of  panels  on  group  J,  non- 
dimensionalised  with  respect  to  the  local  chord,  for  thesidecontaining 
nodes  JINODE(J)  and  J3NODE(J).  TCORDl(J,l)  =  0  0  at  the 
leading  edge  of  the  group  and  TCORDl(J,NC(J)+ 1)  =  1  0  at 
the  trailing  edge.  This  line  of  input  data  is  omitted  if  NCEVEN  =  1 . 

TCORD2(J,K)  ordinates  of  the  chordwise  edges  of  panels  on  group  J,  non- 
dimensionalised  with  respect  to  the  local  chord,  for  the  side  con¬ 
taining  nodes  J2NODE(J)  and  J4NODE(J).  TCORD2(J,l)  =  0  0 
at  the  leading  edge  of  the  group  and  TCORD2(J,NC(J)-(- 1)  =  10 
at  the  trailing  edge.  This  line  of  input  data  is  omitted  if 
NCEVEN  =  I  or  ITCORD  =  0. 

TSORDS(J,K)  ordinates  of  the  spanwise  edges  of  columns  of  panels  on  group  J, 
nondimensionalised  with  respect  to  the  span  of  group  J. 
TSORDS(J.I)  =  0  0  at  the  streamwise  side  containing  node 
JINODE(J),  and  TSORDS(J.NS(J)-f  1)  =  1  0  at  the  streamwise 
side  containing  node  J2NODE(J).  This  line  of  input  data  is 
omitted  if  NSEVEN  =  I. 

NGRPSA  number  of  source  groups,  together  with  any  defined  images,  to  be 

included  in  the  analysis.  The  groups  for  analysis  are  chosen  from 
the  set  of  groups  that  have  been  defined  earlier  in  the  data  file. 

1  ^  NGRPSA  ^  50. 

IGRPSA(K)  list  of  NGRPSA  group  numbers  identifying  the  subset  of  groups 

to  be  included  in  the  analysis  from  the  set  of  groups  defined 
previously.  The  numbers  may  be  in  any  sequence. 

N  MODE  A  number  of  modes  which  are  to  be  used  in  the  analysis. 

I  ^  N  MODE  A  20. 

IMODl :A(K|  list  of  NMODEA  mode  numbers  that  identifies  the  subset  of 
modes  for  which  the  analysis  is  to  be  carried  out.  The  numbers 
m  the  list  refer  to  the  mode  numbers  defined  in  SUBROUTINE 
MODES.  The  numbers  may  be  in  any  order. 

I  he  d.it.i  tile  IS  prim.irils  concerned  with  information  related  to  the  geometry  of  the  con- 
tieur.iiii'ii  heme  .m.il\ee».l  In  order  to  obtain  pressures  and  generalised  forces  it  is  necessary  to 
MippU  ,1  Si  HKOl  IIM  MODES  which  defines  the  modal  data  to  be  used  by  program 
MKIitRII’  I  he  mod.il  data  consls;^  (>f  nondimensional  displacements  and  slopes  at  the 
I  . 'I  m.ilvi  ,.'11,1^,1(11111  pointe  anil,  in  order  to  calculate  the  generalised  forces,  the  non- 
,iinieif.j,  ii,il  ,topl,uelllellI^  .it  the  lift  pi'ints  must  also  be  defined.  This  information  must  be 
o.'o.i  n  ,n.  lle.i  (  (  >\1M<)\  MODAIA  asfollows: 

\K  \  \1I  I  I  K  (NP  :m).  NM  20) 

(  (i\l\ln\  MODAI  A  NOMODE.  D71.(NP,NM).  DZD(NP,NM).  DAD(NP.NM) 
I  '  111  ih'i  I--.',  i.iteO  wiili  (he  ,ibo\e  COMMON  block  are  defined  as  follows: 

Nl’  m.iMiiuim  number  of  panels  which  can  be  used.  This  is  currently 

'•et  to  NP  2W). 

NM  m.isiniuni  luiniber  of  modes  for  which  pressures  can  be  calculated. 

I  hl^  le  currentK  set  to  NM  20. 


IGRPSA(K) 


NMODEA 


IMODEAlK) 


NOMODI 


loi.il  number  of  moiles  ih.it  are  defined  where  NOMODE  C  NM. 
the  inpiii  list  of  nioiles  lo  he  used  in  the  analysis  defines  which 
modes  .ire  to  be  used 


DZL(N,M)  nondimensional  displacement  at  the  lift  point  on  panel  N  in 

mode  M. 

DZD(N,M)  nondimensional  displacement  at  the  normalwash  collocation 

point  on  panel  N  in  mode  M. 

DAD(N,M)  slope  at  the  normalwash  collocation  point  on  panel  N  in  mode  M. 

The  nondimensional  displacements  are  normalised  with  respect  to  the  reference  length.  Panel  N 
corresponds  to  the  panel  in  the  /’th  row  and  y’th  column  of  the  /k’th  group.  N  is  given  by  : 

N  -  INDEXS(K)  +  (J  1)*NC(K)  +  1  1  ^  J  ^  NS(K)  H.l 

I  ^  J  ^  NC(K) 

and  the  vectors  INDEXS,  NC  and  NS  are  described  in  the  following  section. 

In  order  to  facilitate  the  preparation  of  the  modal  data,  the  specifications  of  a  number  of 
labelled  COMMON  blocks  are  given  below.  These  COMMON  blocks  contain  data  related  to 
the  geometry  of  the  configuration,  as  well  as  information  associated  with  the  organisation  of 
this  data,  and  are  given  below  (note  that  NP  was  defined  previously): 

PARAMETER  (NE  50) 

COMMON,  AERPAR  GNU.  AM,  RL,  AM2,  BETA2.  BETA.  RL2 
COMMON  XYZLFT  XL(NP.4),  YL(NP,4),  ZL(NP,4) 

COMMON  XYZDSH  XD(NP.4),  YD(NP,4),  ZD(NP,4) 

COMMON,  INDEX  INDEXS(NE),  INDEXA(NE) 

COMMON  ELDATA,NOELEM,  NC(NE),  NS(NE),  IREF(NE,4).  lELUSE(NE) 
COMMON, DIHEDL  D1HEDL(NE,4) 

The  definitions  of  the  variables  are  : 


NE 

maximum  number  of  groups  which  can  be  defined. 

This  is  currently  set  to  NE  =  50. 

GNU 

frequency  parameter 

AM 

Mach  number 

RL 

reference  length 

AM2 

AM  »  AM 

BETA2 

1  0  AM  »  AM 

BETA 

SQRT  (BETA) 

RL2 

RL  *  RL 

XL(N.IR), 

YL(N,IR). 

ZL(N.IR) 

dimensional  (X.  Y.  Z)  coordinates  of  the  lift  point  on  the  /V'th 
panel  m  the  //?'th  reflection  state.  /R  1  refers  to  the  source 

group  (see  Section  8.3). 

XD(N.IR). 

YD(N,IR). 

ZD(N.IR) 

dimensional  (X.  Y,  Z)  coordinates  of  the  normalwash  collocation 
point  on  the  /V'th  panel  in  the  IR'lU  reflection  state. 

INDEXS(K) 


pointer  indexing  the  position  of  the  first  panel  associated  with 
the  A'th  group.  This  pointer  refers  to  all  source  groups. 


INDEXA(K) 


pointer  indexing  the  position  of  the  first  panel  associated  with  the 
/f’th  group  used  in  the  analysis. 

NOELEM  number  of  groups  (elements)  defined  in  the  data  file. 


NC(K) 

NS(K) 

1REF(K..1R) 


lELUSE(K) 


number  of  chordwise  panels  for  the  Af'th  group. 

number  of  spanwise  panels  for  the  A"th  group. 

number  defining  the  sense  of  positive  pressure  on  the  A^th  source 
group  (IR  ^  I)  and  any  image  groups  (IR  =  2,'  3  or  4).  See  the 
definition  of  IREF  used  in  relation  to  the  input  data  file. 

this  variable  indicates  whether  a  particular  group  (clement) 
number  has  been  used.  If  1ELUSE(K.)  =  0  then  the  A'th  group 
(element)  has  not  been  defined.  If  lEHJSE(K)  =  1  then  the 
A’th  group  has  been  defined. 


D1FIEDL(K,1R)  this  is  the  dihedral  angle  (in  radians)  of  the  A'th  group.  IR  =  1 
refers  to  the  source  group  and  IR  —  2,  3  and  4  refer  to  any  image 
groups  that  were  created. 


Let  us  consider  the  iV'th  panel  oscillating  in  mode  M.  The  coordinates  of  the  normalwash 
collocation  point  may  be  given  as  (.v.s,  vn,  cn),  and  by  applying  Equation  6.1  the  displacement 
normal  to  the  surface  at  the  collocation  point  may  be  expressed  in  the  following  form  : 


•?m(-v.n,.i'x,  cn.  /)  =  //m(.vn,,i'.n,  rN)t’'‘'“ 


H.2 


The  coordinates  of  the  normalwash  collocation  point  on  the  A'"th  source  panel  are  given  by  the 
following  equation  : 


(.VN. .I  N.  rs)  (T/)(,V.|).  )72(.\'.I),ZZ)(AJ))  H.3 

Using  Equation  6.3  the  nondimensionui  normalwash  may  be  expressed  as 


X/  >.(  /'tl'A-  'N-  -n)  1  + 'k7M(.V.\.  )■ 

Nx  -n) 

H.4 

l)\  j 

DAD(N.M)  f /A  .  DZD(N.M) 

H.5 

In  order  to  obtain  the  generalised  forces  it  is  necessary  to  calculate  the  nondimensional 
displacement  at  the  mid-span  quarter-chord  point  of  the  A"th  source  panel.  (.v,\.  vn.  cn),  in  the 
iV/  th  mode.  The  nondimensional  displacement  may  be  determined  from  the  following  equation  : 


./m(  V M .  r M.  CM )  /m(  ,\7,(  A.  I ).  >  /,(  ,V.  I ).  ZL(  ,V.  I ) ) 


H.6 


DZUNM)  H.7 

Although  it  is  customary  to  defme  OZL  in  the  manner  depicted  above,  this  convention  need 
not  be  loDovved  il  lorces  other  than  generalised  forces  are  required. 

A  sign  convention  exists  for  the  specification  of  the  modal  displacements  and  slopes  used 
in  the  calculation  ol  pressures  and  generalised  forces.  The  vlisplacements  perpendicular  to  the 
surface.  l)/l)  and  0/1..  are  taken  to  be  positive  when  Ihey  are  in  the  same  direction  on  a  group 
as  the  direction  ol  positive  pressure  fi'r  that  group.  The  slope  at  the  normalwash  colloeation 
point.  DAI),  is  positive  if  tlie  displacement  at  that  point.  D/D.  is  increasing  in  the  streamvvisc 
( A  t  direction.  An  example  illustrating  the  application  of  the  above  convention  is  given  in  Appen¬ 
dix  I. 


I 

I 


APPENDIX  1 


Example  of  input  preparation  for  the  analysis  of  Stark's  T-tail 

As  an  illustrative  example,  the  preparation  of  the  input  data  file  and  SUBROUTINE 
MODES  for  Stark’s  T-tail  will  be  described.  The  configuration  is  assumed  to  be  oscillating  with 
a  frequency  parameter  of  0  -6,  based  on  a  reference  length  of  one  unit,  in  an  airstream  of  Mach 
number  0-8. 

The  data  file  for  the  panel  distribution  shown  in  Figure  17  is  given  in  Table  13.  A  total  of 
7  nodes  were  used  to  define  the  configuration.  The  positions  of  the  nodes  can  be  deduced  by 
using  the  node  coordinates  in  conjunction  with  Figure  16.  Two  groups  were  defined,  with 
group  I  corresponding  to  the  starboard  half  of  the  stabiliser  and  group  2  corresponding  to  the  fin. 

Due  to  symmetry,  only  the  starboard  half  of  the  stabiliser  is  modelled.  Nodes  1  and  2  define 
the  leading  edge  and  nodes  3  and  4  define  the  trailing  edge  of  the  stabiliser.  IREF(1,1)  =  I 
indicates  that  the  direction  of  positive  pressure  on  the  source  group  is  in  the  negative  r-direction. 
This  agrees  with  the  topology  of  the  group  (see  Section  8.2).  The  effect  of  the  port  half  of  the 
stabiliser  is  included  by  setting  IREF(I,2)  -  1.  This  creates  an  image  of  the  source  group  by 
reflection  in  the  .v-c  plane. 

The  positive  sense  of  pressure  on  the  image  is  in  the  positive  r-direction,  which  once  again 
corresponds  with  the  topology  of  the  group.  The  combination  of  positive  pressures  on  the  port 
and  starboard  halves  of  the  stabiliser  corresponds  to  antisymmetric  motion  of  the  stabiliser. 

The  stabiliser  is  divided  into  1 1  panels  chordwise  and  10  panels  spanwise,  with  the  chord- 
wise  panel  distribution  being  the  same  on  both  sides  of  the  half  stabiliser.  The  chordwise  and 
spanwise  ordinates  of  the  panel  sides  are  listed  in  lines  16  and  17  of  the  data  file. 

The  leading  edge  of  the  fin  is  defined  by  nodes  2  and  5  and  the  trailing  edge  by  nodes  6  and  7. 
Since  IREF(2,I)  -  —I,  the  sense  of  positive  pressure  on  the  fin  is  in  the  negative  (-direction 
which  is  opposite  to  that  defined  by  the  topology  of  this  source  group.  Note  that  no  image 
surfaces  are  created.  The  fin  is  divided  into  9  panels  chordwise  and  10  panels  spanwise,  and  the 
nondimensional  ordinates  of  the  panel  sides  are  listed  in  lines  22  and  23  of  the  data  file. 

For  this  analysis,  both  the  stabiliser  and  fin  are  included,  which  is  indicated  in  lines  24  and  25. 
Three  modes  are  to  be  used  in  the  calculation  of  pressures  and  generalised  forces.  The  numbers 
1.  2  and  3  in  the  last  line  refer  to  the  modes  defined  in  SUBROUTINE  MODES. 

Since  IPLOT  -  0  (see  line  2)  a  PAFEC  data  file  for  use  in  PIGS  will  not  be  created.  By 
setting  IPRNTP  I  the  print  out  of  pressures  associated  with  the  three  modes  will  be  presented 
in  a  columnar  format. 

A  SUBROUTINE  MODES  is  required  to  generate  the  modal  data  and  a  detailed  description 
of  the  format  of  this  subroutine  is  given  in  Appendix  H.  A  listing  of  the  FORTRAN  source 
code  for  the  subroutine  used  to  calculate  modal  displacements  for  Stark's  T-tail  is  presented  in 
Table  14. 

As  described  in  Section  9.2.  the  T-tail  is  assumed  to  be  oscillating  in  three  rigid  body  modes. 
These  are  yawing  about  a  vertical  axis  through  the  centre  of  the  root  chord  of  the  fin  (positive 
nose  right),  sidesway  (positive  left)  and  rolling  about  the  fin-stabiliser  intersection  (positive  right 
stabiliser  down).  These  modes  are  expressed  analytically  in  Equations  9.5  to  9.10,  and  when 
they  are  used  in  Equations  H.4  to  H.7  the  following  equations  are  obtained  for  the  stabiliser ; 


DZD(N.I) 

DAD(N.I) 

DZl.(N.I)  00 

l.l 

DZD(N.2) 

DAD(N.2) 

DZI.(N,2)  0  0 

1.2 

D/'D(N.3)  MJIN.I) 


1.3 


DAD(N,3)  =  0  0  1.4 

DZL(N,3)  =  YL(N,1)  1.5 

where 

N  =  INDEXS(K)  +  (J-1)*NC(K)  +  I  1.6 

1  <  I  <  NC(K)  1.7 

1  <  J  ^  NS(K)  1.8 

and  K  =  1  for  the  stabiliser  group.  The  modal  data  for  the  fin  is  as  follows : 

DZD(N,1)  =  3  0  *(XD(N.l)  +  0- 15577)  1.9 

DAD(N,1)  =  30  1.10 

DZL(N,1)  =  3  0  *(XL(N,l)  +  0- 15577)  1.11 

DZD(N,2)  =  1  0  1.12 

DAD(N,2)  =  00  1.13 

DZL(N,2)=10  1.14 

DZD(N,3)  =  -ZD(N,1)  1.15 

DAD(N,3)  =  00  1. 16 

DZL(N,3)  -  -ZL(N,1)  1.17 

where  N  is  obtained  from  Equations  1.6  to  1.8  by  substituting  K  =  2  for  the  fin  group. 


Note  that  the  displacements  DZL  and  DZD  are  positive  in  the  same  direction  as  the  sense 
of  positive  pressure  that  is  defined  for  the  panels  in  a  given  group.  The  slope  DAD  is  positive 
if  the  displacement  is  increasing  in  the  streamwise  .v-direction. 

Figure  20  illustrates  the  mode  shapes  for  the  three  modes  used  for  Stark’s  T-tail.  The  4  sign 
indicates  that  the  displacement  is  in  the  same  direction  as  the  sense  of  positive  pressure  acting  on 
the  group,  and  the  sign  indicates  that  the  displacement  is  in  the  direction  opposite  to  that  of 
the  positive  pressure.  Hence,  where  the  +  sign  applies  the  displacements  are  positive  quantities 
and  where  the  sign  applies  the  displacements  are  negative  quantities.  It  is  seen  that  when  the 
analytic  expressions  for  the  mode  shapes  are  combined  with  the  topology  of  the  groups  (as  de¬ 
fined  in  the  input  data  file)  then  the  displacements  depicted  in  Figure  20  are  obtained. 


TABLE  1 


Values  of  Generalised  Airforces,  Qp,„  for  the  AGARD  Wing-Tailplane  Configuration  (for  h  0) 

as  Obtained  by  Different  Workers 


00 

o  o 

II  11 

Qw 

Davies 

Albano 

Perkinson 

Rodden 

1 

Albano 

Rodden 

Geising 

Kalman 

Rodden 

Farrell 

Present 

Method 

Qii 

0-4403 

0-4425 

0-4554 

0-4401 

0-4377 

0-4557 

/IQii 

359-9 

359-9 

359-9 

359-9 

360-0 

360-0 

1012  1 

0-6202 

0-6121 

0-6655 

0-6557 

0-6457 

0-6652 

A:  —  0 

.^012 

180-0 

180-0 

1 

180-0 

180-0 

180-0 

180-0 

1  021  1 

0-1046 

0-1054 

^  0  1107 

0- 1044 

0- 1049 

0-1083 

Z021 

180-3 

180-3 

180-3 

180-3 

180-0 

180-0 

022 

0-1759 

0-1954 

0-2237 

0-2126 

0-2184 

0-2261 

/  ^22 

180-2 

180-2 

180-2 

1 

180-2 

180-0 

180-0 

011 

1  -  5865 

1 -6022 

1 -5496 

1 -5688 

1  -4212 

1-5713 

..  011 

314-2 

314-4 

;  311  2 

i 

310-7 

312-3 

310-8 

012 

0-9180 

0-8910 

0-9081 

0-9495 

0-8752 

0-9482 

k  ~  1-5 

/012 

265-5 

266  3 

267-2 

265-4 

265-0 

265-8 

1021 

1 -0043 

1 -0099 

1 -0550 

1-0511 

0-9507 

1 -0653 

./021 

291-7 

291  -4 

287-2 

288-5 

289-9 

289-0 

022 

1 -2845 

1 -2386 

1 -2144 

1-2719 

1-1448 

1 -2745 

/  022 

294-7 

294-3 

293-7 

292  9 

294-1 

293-2 

TABLE  2 


Values  of  Generalised  Airforces,  Qp,|,  for  the  AGARD  Wing-Tailplane  Configurations  as  Obtained 

by  Different  Workers 


h  =  0'6 

M  =  0-8 

0p<i 

Davies 

Albano 

Perkinson 

Rodden 

Farrell 

Present 

Method 

11 

O 

1  011  ! 
Z011 

0  1470 
359-8 

0-1490 

359-8 

0-1374 

360-0 

0-1432 

360-0 

1  012  i 
Z012 

0-6402 

180-1 

0-6312 

180-1 

0-6661 

180-0 

0-6868 

180-0 

021  1 

/:02i 

0-2404 

180-1 

0-2405 

180-1 

0-2527 

180-0 

0-2669 

180-0 

022  i 
./  022 

0-1619 

180-3 

0- 1817 
180-2 

0-1958 

180-0 

0-2117 

180-0 

A-  -  1  -5 

i  011  i 
Z011 

1-1009 

301  3 

1-1200 

301  -3 

0-9780 

299-9 

1 -0837 
298-3 

'  012  1 
Z012 

1-1342 

251  -7 

1-1128 

251  -9 

1 -0735 
249-5 

1-1583 

250-0 

021  1 
/021 

0-9072 

278-0 

0-9122 

277-7 

0-8259 

276-9 

0-9239 

275-8 

1 -3867 
289-2 

1 -3397 
288-6 

1 -2082 
287-2 

1 -3541 
286-4 

TABLE  3 


Generalised  Airforces  for  the  AGARD  Wing-Tailplane  Configuration  (with  h  0-6,  M  ^  0-8) 

for  a  Range  of  Frequency  Parameters 


TABLE  4 


Generalised  Airforces  for  the  AGARD  Wing-Tail  Configuration  for  Various  Values  of  Vertical 

Separation,  h,  Between  Wing  and  Tail 


II  II 

o  — 

00  C/i 

1011  1 

LQn 

:  Qiz  1 
/_Qi2 

1  Qn  i 

IQii 

Qii  ^ 

LQ-ii 

h 

Davies 

Present 

Method 

Davies 

Present 

Method 

Davies 

Present 

Method 

Davies 

Present 

Method 

0 

1-5865 

314-2 

1  -5713 
310-8 

0-9180 

265-5 

0-9482 

265-8 

1  -0043 
291  -7 

1  -0653 
289-0 

1-2845 

294-7 

001 

1-5519 

313-3 

1-5514 

310-7 

0-9332 

264-0 

0-9520 

264-9 

0-9973 
290  4 

1  -0574 
288-7 

1  -2890 
294-0 

1  -2734 
292-8 

004 

1  -4773 

31 1  -4 

1  -4777 
309-3 

0-9640 

261-2 

0-9789 
261  -7 

0-9843 

m-1 

1  0347 
286-6 

1  -  2992 
292-6 

1  -2782 
291  -3 

01 

1  -3730 
308-7 

1-3605 

306-1 

1-0040 

258-0 

1-0303 

257-5 

0-9662 

284-3 

1  -0042 
282-8 

1  -3151 
291  -3 

1  -2946 
289-2 

0-2 

1  -2624 
305-9 

1  -2461 
303-0 

1  -0479 
255-2 

1  -0779 
254-1 

0  9428 
281-4 

0-9724 

279-5 

1  -3376 
290-4 

1  -3138 
287-8 

0-3 

1-1948 

304-0 

1-1787 
301  -2 

1  -0790 
253-7 

1  - 1072 
252-3 

0-9271 

279-8 

0-9520 

277-8 

1-3559 

290-0 

1  -3270 
287-2 

0-4 

1-1511 

302-8 

t-1349 

299-9 

1  - 1024 
265-8 

1-1284 
251  -2 

0-9174 

279-0 

0-9383 

276-8 

1-3695 

289-7 

1-3376 

286-9 

05 

1  1216 
301  -9 

1-1048 

299-0 

1-1204 

252-2 

1-1450 

250-5 

0  9113 
278-4 

0-9295 

276-2 

1  -3795 
289-5 

1-3468 

286-6 

0-6 

1-1009 
301  -3 

1  -0837 
298-3 

I -1342 
251  -7 

1  - 1 583 
250-0 

■11 

B 

0  9239 
275-8 

1  -3867 
289-2 

1  -3541 
286-4 

TABLE  6 


Generalised  Airforces  for  Stark’s  T-tail  as  Obtained  by  Different  Workers 


k  = 
M  = 

0-6 

0-8 

k  =  0-9 
M  -  0-8 

Stark 

Davies 

Zwaan 

Kalman 

Rodden 

Geising 

Farrell 

Present 

Method 

Stark 

Farrell 

Present 

Method 

3  0826 
258-8 

3-2421 

260-6 

3  -  2873 
259-8 

3-3527 
261  0 

2 • 7470 

260  6 

3-0965 

260-5 

4-8800 

264-9 

4  - 1 775 
266-2 

4-8056 

265-5 

0  3399 
328-5 

0  3475 
327-4 

0  3431 
329-4 

0-2799 
328  1 

0-3214 
328  - 1 

0  7020 
330  6 

0-6016 
333  -  7 

0-7042 

332-7 

0-1695 

62-4 

0-1859 

61-5 

0- 1865 
60-8 

0-1830 
61  9 

0- 1680 

61  6 

0-1828 

60  9 

0  3214 
52  0 

0-2971 
50  4 

0-3358 

49-7 

4-4628 

210-7 

4-5670 

211-1 

4  5650 
210  9 

4  6612 
212-3 

4-4291 

211-0 

4  6085 
21  1  0 

5-3736 

220-8 

5-1 126 
220- 1 

5-4472 

221-2 

0-7758 
281  -2 

0-7930 
282- 1 

0  7936 
281  -9 

0-8134 

283-3 

0-7815 

282-3 

0-8072 

282-2 

1  2419 
286-5 

1  -2185 

288-3 

1  -2822 
287-6 

0-2103 

298-2 

0  2191 
299-0 

0-2183 

298-9 

0-2269 

298-3 

0-2282 

296-6 

0-2330 

297-8 

0  3692 
310-2 

0-3882 

307-6 

0-4055 
309  0 

1-0772 

222-7 

1-1154 

224-7 

1-1006 

224-6 

1-1383 

224-9 

1-1483 

225-3 

1-1686 

224-9 

1 -3760 
232- 1 

1  4446 
236  6 

1  -4903 
235-3 

0-1810 

297-7 

0- 1897 
299-8 

0-1874 

299-7 

0-1948 

299-8 

0-2000 

300-1 

0-2022 

299-7 

0-3214 

305-8 

0-3484 

308-5 

0-3557 

307-6 

0-3319 

289-8 

0-3349 

289-6 

0-3349 

289-6 

0-3523 

289-4 

0-3590 

288-8 

0-3617 

289-4 

0-5484 

298-1 

0-5794 

296-8 

0-5910 

297-5 

TABLE  7 


Comparison  of  Results  for  Stark’s  T-tail  Obtained  with  Different  Panel  Distributions.  The  Top 
Table  Indicates  the  Panel  Distributions  Used 


TABLE  8 


Generalised  Airforces  for  the  ONERA  Horizontal  Wing  and  Tail  Configuration  for  H  1/8 

and  A  2 


W=0-30 

A:=0-3856 

M=Q 

k=Q 

-45 

•2436 

A/=0-65 

A:=0-1513 

A/=0-80 

A-  =0-1112 

Davies 

Present 

Method 

Davies 

Present 

Method 

Davies 

Present 

Method 

Davies 

i 

^  Present 
Method 

14-209 

349-9 

14-068 

354-8 

14-786 

353-2 

14-850 

355-6 

16-248 

354-0 

16-396 

355-7 

18-346 

353-4 

18-620 

354-9 

0-214 

311-7 

0-244 

307-6 

0-234 

311  3 

0-260 

309-6 

0-220 

301  -2 

0-241 

300-9 

. 

0-181 

273-9 

0- 196 
273-9 

8-057 

1 16-4 

8-252 

1 12-5 

8  518 
138-4 

8  675 

136  0 

9-897 

152-3 

10-039 

150-8 

12-066 

157-4 

12-293 

156-3 

14-130 

351-5 

14-041 

355-1 

14-743 

353-4 

14-804 

355-8 

16-198 

354-2 

16-342 

355-9 

18-301 

353-6 

18-572 

355-1 

1  066 

147-7 

1  -221 
143-5 

1  150 
144-9 

1  -290 
143-0 

1  -081 
135-8 

1  194 
135-3 

1 

0-872 

1 14-7 

0  957 

1 14-5 

29  102 
296  3 

30  196 
292-0 

H 

31-464 
316  1 

34  •  703 
333  •  3 

35-716 
331  -6 

41 -054 
339-2 

42-315 

337-9 

TABLE  9 


Generalised  Airforces  for  the  ONERA  Horizontal  Wing  and  Tailplane  Configuration  for  H  1/8 

and  A  1/2 


//=l/8 

A=l/2 

Af=0-30 

A-=0-3856 

Af=0-45 

A- =  0-2436 

A/-0-65 

^=0-1513 

M=0-80 

k=-Q\\\l 

^ii'i 

Davies 

Present 

Method 

Davies 

Present 

Method 

Davies 

Present 

Method 

Davies 

Present 

Method 

033 

13-847 

13-797 

14-417 

14-503 

15-889 

16-058 

18-075 

18-381 

033 

352-7 

356-39 

354-3 

356-8 

354  9 

356-7 

354-4 

355-8 

03  4 

0-900 

0-893 

0-928 

0-929 

0-857 

0-852 

0-702 

0-693 

:!  03  4 

324-2 

321 -89 

326  4 

3'5-2 

318  3 

320-6 

305  - 1 

302-7 

043 

7-934 

8-118 

8-415 

8-561 

9-732 

9-865 

1 1  -813 

12-025 

043 

156- 1 

152-1 

164- 1 

161-7 

168  9 

167-3 

170- 1 

169-0 

0.4, 

13-279 

13-163 

1 3  -  794 

13-812 

15-200 

15-298 

17  345 

17-562 

0.4 

353-8 

357-7 

355  1 

357  8 

355-6 

357-4 

355-1 

356-7 

0.4 

6-381 

6-791 

6-641 

7  106 

6  446 

6-886 

5  -  666 

6-067 

014 

165-5 

164  0 

164-5 

163  8 

161  -2 

161  -  1 

153-4 

153-7 

033 

28-082 

29-073 

29-469 

30-370 

33-518 

34-429 

39-664 

40-805 

033 

334-9 

330-7 

343-7 

341-1 

349  3 

347-6 

351  -4 

350-1 

/ 

/'i 

Ci 

1 

0  984807753 

0- 007642844 

■) 

0  860025404 

0  469841270 

3 

0  (.42787610 

0- 003 167376 

4 

0  342020143 

5-610810220 

5 

0  0 

14-139682541 

(1 

0  342020143 

5-610810220 

7 

0  (.42787610 

0-003167376 

0  866025404 

0-469841270 

9 

0  984807753 

0-007642844 

TABLE  13 


Listing  of  data  file  used  for  the  analysis  of  Stark’s  T-tail  configuration 


0-6  0-8  10 
0  1 
7 

1  0-341  10  0-0 

2  0  000  0  0  0  0 

3  0-813  1-0  0-0 

4  0-938  0-0  0-0 

5  0-801  0-0  1-0 

6  0-820  0-0  0-0 

7  0-489  0-0  1-0 

2 

STABILISER 

I  12  3  4 
110  0 

II  10  0  0  0 

0-0  -05245  13113  -21855  -36716  -50704  -65565  -74307  82175  -87420  -941  1-0 


0-0  -04  -10 

FIN 

2  2  5  6  7 

- 18 

-30 

50 

•70 

-82 

-90 

•96 

1  -0 

10  0  0 

9  10  0  0  0 

00  -06  15 

-25 

-42 

-58 

•75 

•85 

•94 

1  -0 

0-0  04  10 

•  18 

-30 

50 

70 

•82 

•90 

•96 

1  -0 

2 

I  2 
3 

1  2  3 


TABLE  14 


Listing  of  FORTRAN  source  code  for  the  subroutine  used  to  generate  the  modal  data  for 

Stark’s  T-tail  configuration 


1  C 

2  ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

3  C 

4  SUBROUTINE  MODES 

5  C 

6  C  THIS  SUBROUTINE  CALCULATES  THE  MODES  FOR  THE  UNSTEADY 

7  C  AERODYNAMICS  PROGRAM. 

8  C 

9  C  THE  MODES  ARE  FOR  STAPK'S  T-'l  AlL. 

10  C 

11  C  MODE  1:  YAW 

12  C  MODE  2:  SIDESLIP 

13  C  MODE  3:  ROLL 

14  C 

15  ccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

16  C 

17  PARAMETER  (NP  260,NE-50,NM  -  20) 

18  C 

19  COMMON  AERPAR/GNU, AM, RL 

20  CO  M  M  O  N ,  X  Y  Z  L  FT  X  L(  N  P.4),  Y  L(  N  P,4).Z L(  N  P.4 ) 

21  COMMON  XYZDSH/XD(NP,4),YD(NP,4),ZD(NP.4) 

22  COMMON  ELDATA  NOELEM,NC(NE),NS(NE),IREF(NE,4) 

23  COMMON  INDEX  / INDEXS(NE).INDEXA(NE) 

24  C 

25  COMMON  MODATA,  NOMODE.DZL(NP.NM),DZD(NP,NM).DAD(NP,NM) 

26  C 

27  C  NUMBER  OF  MODES  DEFINED 

28  C 

29  NOMODE -3 

30  C’ 

31  C  MODES  FOR  STABILISER:  GROUP  I 

32  C 

33  DO  10  N  INDEXS(I)4- l,INDEXS(l)  +  NC(l)*NS(l) 

34  C 


35 

36 

37 


38 

39 

40 

41 

42 

43 

44 

45 

46 

47 


C 


C 


( 


DZD(N.I)  00 
DAD(N,I)  00 
DZLIN.I)  00 

DZD(N.2)  0  0 
DAD(N.2)  00 
DZL(N.2)  0  0 

DZD(N,3)  ^YD(N.I) 
DAD(N.3)  0  0 
DZL(N.3)  -YL(N,I) 


10  CONTINUE 


48  C 

49  (  MODES  FOR  FIN:  GROUP  2 

50  C 


TABLE  14  (Continued) 


51  DO  20  N  =  INDEXS(2)+1,INDEXS(2)  +  NC(2)*NS(2) 

52  C 

53  DZD(N,  1 )  =  +  3  •0*(XD(N.  1 )  +  0  ■  1 5577) 

54  DAD(N,l)=+30 

55  DZL(N,1)  =+3-0*(XL(N.l)  +  0- 15577) 

56  C 

57  DZD(N,2)  =  +  10 

58  DAD(N,2)=0  0 

59  DZL(N,2)-  +  10 

60  C 

61  DZD(N,3)=-ZD(N,1) 

62  DAD(N,3)=00 

63  DZL(N,3j=  ZL(N.I) 

64  C 

65  20  CONTINUE 

66  C 

67  RETURN 

68  END 


FIG.  1:  EXAMPLE  OF  SURFACE  REPRESENTATION  BY  PANELS  AND  THE  FEATURES 
OF  THE  PANEL  GEOMETRY. 


Lift  point 


Sending  panel  s 


FIG.  2 


B 


Control  point 


Point  A:  (X,  Y,Z) 
Point  B;  (X,Y,  Z) 
Point  C:  (X,  Y,Z) 


(xa,  Va, Za) 
(xb,  Vb,  Zb) 
(x,  y,  z) 


Panel  dihedral  angles  :  and  y^ 


DIAGRAM  SHOWING  A  GENERAL  COMBINATION  OF  A  SENDING  AND 
RECEIVING  PANEL, 


COMPARISON  OF  NONPLANAR  AND  PLANAR  PRODUCT  OF  THE  STEADY 
KERNEL  FUNCTION  AND  r2,  Kor2,  for  M  =  0.6. 


Parameters  used  in  expression  for  Kor2 : 
=  60°  7j  =  -30° 

M  =  0.6 

xi  =  1.0 


Vi  =  COS7j  t 
z,  =  sin7jt  +  0.04 


FIG.  5:  VARIATION  OF  Kor?  FOR  A  SENDING  AND  RECEIVING  PANEL  COMBINATION 
WHERE  THE  PANELS  ARE  PERPENDICULAR  TO  EACH  OTHER. 


FIG.  8:  DIAGRAM  SHOWING  NODE  NUMBERING  SCHEME  AND  OTHER  FEATURES 
PERTAINING  TO  THE  DEFINITION  OF  A  GROUP.  THE  NUMBERS  1,  2,  3  AND 
4  REFER  TO  THE  NODES. 


J'th  column 


Node  1  :  (tc,  tg)  =  (0,  0) 
Node  2  :  (t^  tg)  =  (0,  1) 
Nodes  :  (tc,  tg)  =  (1,0) 
Node  4  :  (tc,  tg)  =  (1,1) 


FIG.  9;  DIAGRAM  DEFINING  NONDIMENSIONAL  COORDINATE  SYSTEM  FOR  A  GROUP. 
ALSO  SHOWN  IS  THE  PANEL  NUMBERING  CONVENTION  FOR  PANELS  WITHIN 
A  GROUP. 


FIG.  13:  TOPOLOGY  OF  IMAGE  GROUPS,  SUPERSCRIPT  (2),  (3)  AND  (4),  WHEN 
GENERATED  FROM  A  SOURCE  GROUP  (1). 


FIG.  14:  IDEALISATION  OF  T-TAIL  WITH  TWO  AXES  OF  SYMMETRY.  ONLY  TWO 
SOURCE  GROUPS  NEED  TO  BE  DEFINED,  THE  OTHER  SURFACES  BEING 
OBTAINED  BY  REFLECTION.  THE  ARROWS  INDICATE  THE  POSITIVE 
SENSES  OF  PRESSURE  ASSUMED  BY  THE  TOPOLOGY  OF  THE  SOURCE 
GROUPS. 


FIG.  15:  AGARD  HORIZONTAL  WING-TAILPLANE  CONFIGURATION. 


FIG.  16:  STARK'S  SWEPT 


Yaw  of  fin 


Sidesway  of  fin 


J 


Roll  of  fin  and 
stabiliser 


FIG.  20:  SKETCH  OF  MODE  SHAPES  USED  IN  THE  ANALYSIS  OF  STARK'S  T-TAIL. 

THE  +  AND  -  SIGNS  INDICATE  THAT  THE  DISPLACEMENTS  ARE  IN  THE  SAME 
OR  OPPOSITE  DIRECTION  TO  THE  SENSE  OF  POSITIVE  PRESSURE  FOR  THE 
GROUP. 
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